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In brief 

When and how echolocation first evolved 
in bats is contentious. Hand et al. report 
on an uncrushed, 50-million-year-old bat 
skull that provides significant new data to 
inform the debate. Measurements of the 
inner ear of the French stem bat Vielasia 
suggest that advanced (laryngeal) 
echolocation may have evolved before 
the crown bat radiation. 
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Highlights 
e Anew bat species is described from a 50-million-year-old 
cave deposit in France 


e Its well-preserved fossils include the oldest uncrushed skull 
of a bat yet known 


e This stem bat appears to have been capable of advanced 
(laryngeal) echolocation 


e The fossils suggest that advanced echolocation predates the 
crown bat radiation 
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SUMMARY 


Bats are among the most recognizable, numerous, and widespread of all mammals. But much of their fossil 
record is missing, and bat origins remain poorly understood, as do the relationships of early to modern bats. 
Here, we describe a new early Eocene bat that helps bridge the gap between archaic stem bats and the hy- 
perdiverse modern bat radiation of more than 1,460 living species. Recovered from ~50 million-year-old cave 
sediments in the Quercy Phosphorites of southwestern France, Vielasia sigei’s remains include a near-com- 
plete, three-dimensionally preserved skull—the oldest uncrushed bat cranium yet found. Phylogenetic ana- 
lyses of a 2,665 craniodental character matrix, with and without 36.8 kb of DNA sequence data, place Vielasia 
outside modern bats, with total evidence tip-dating placing it sister to the crown clade. Vielasia retains the 
archaic dentition and skeletal features typical of early Eocene bats, but its inner ear shows specializations 
found in modern echolocating bats. These features, which include a petrosal only loosely attached to the ba- 
sicranium, an expanded cochlea representing ~25% basicranial width, and a long basilar membrane, collec- 
tively suggest that the kind of laryngeal echolocation used by most modern bats predates the crown radia- 
tion. At least 23 individuals of V. sigei are preserved together in a limestone cave deposit, indicating that cave 
roosting behavior had evolved in bats by the end of the early Eocene; this period saw the beginning of sig- 
nificant global climate cooling that may have been an evolutionary driver for bats to first congregate in caves. 


INTRODUCTION 


Bats are one of the most speciose and ecologically diverse 
mammalian orders, representing one-fifth of all mammals alive 
today (>1,460 species'). They are the only mammals to use pow- 
ered flight, and many bats (but not flying-foxes or Old World fruit 
bats) use echolocation calls produced in the larynx to orientate in 
darkness and find food. These two key biological innovations 
make bats morphologically highly distinctive in their skeletons 
and skulls, and molecular data suggest they diverged from other 
mammals more than 60 million years ago (mya).” © 

When and how flight and echolocation evolved in bats is 
contentious. It has been argued that powered flight evolved 
before echolocation in bats,” ° that echolocation came first, '° 
or that powered flight and echolocation evolved in tandem. '' 
Today, all bats fly, but not all bats echolocate. Old World fruit 
bats and flying-foxes (Pteropodidae, suborder Yinpterochirop- 
tera) do not use laryngeal echolocation, unlike the other ~86% 
of extant bats, which utilize high-frequency sounds produced 
by the larynx to echolocate. There is debate about when and 


how laryngeal echolocation evolved in bats, with two specific hy- 
potheses debated: (1) that laryngeal echolocation evolved inde- 
pendently within the two bat suborders Yinpterochiroptera and 
Yangochiroptera (e.g., Teeling et al.'* and Nojiri et al.) or (2) 
that laryngeal echolocation evolved once in the common 
ancestor of extant bats but was subsequently lost in some yinp- 
terochiropterans (i.e., pteropodids; e.g., Simmons et al.,° 
Springer et al.,'“ Teeling, ® and Liu et al.'°). A range of morpho- 
logical, genomic, anatomical, phylogenetic, developmental, 
physiological, and behavioral data has been used to support 
both alternatives (e.g., Nojiri et al.,'° Liu et al.,'° Eick et al.,'’ Ve- 
selka et al.,'° Davies et al.,'° Teeling et al.,°° Wang et al.,™ and 
Thiagavel et al.”°). 

The bat fossil record also provides significant evolutionary in- 
sights, although it is one of the poorest of any mammalian order, 
with ~80% of the record estimated to be missing.°?*:2* The 
world’s oldest known bat fossil consists of a single, 56 mya lower 
molar from Portugal.°°-*’ Another fifty other early to middle 
Eocene bats are known globally.” These are referred to several 
extinct families, but most are known from fragmentary dental 
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remains, and their relationships to modern bats remain poorly 
understood. Current evidence indicates that these early bats 
fall outside the crown clade,*° °° with the first definitive repre- 
sentatives of modern bat families appearing in the fossil record 
from ~50 mya (MP10).°°34°° 

A few early and middle Eocene stem bats are represented by 
whole-bodied remains preserved in oil shales from lacustrine 
Konservat-Lagerstatten such as Green River, Wyoming, USA, 
and Grube Messel, Darmstadt, Germany. Although flattened, 
they provide fine details of the cranium, postcranium, and soft 
tissues.°°°“° For these rare early bats, a range of echolocation 
capabilities has been argued on the basis of cochlea size, the 
shape of the orbicular process of the malleus, and stylohyal 
morphology™'® 07 —from the absence of laryngeal echoloca- 
tion in the 52 mya Onychonycteris finneyi’'*® (but see Veselka 
et al.'°) to near-modern capabilities in 48 mya Hassianycteris 
and Palaeochiropteryx species.*’ However, it has also been sug- 
gested, on the basis of inner ear embryology in extant bats, that 
stem bats may have used a primitive type of echolocation 
different from the sophisticated system used by modern echolo- 
cating bats.'° For early Eocene bat fossils, common preserva- 
tion artifacts, such as flattening, fracturing, disarticulation, and 
deformation, complicate interpretations of the spatial relation- 
ships of delicate bones and, in turn, echolocation abilities. '°7° 

Here we describe a new, three-dimensionally preserved fossil 
bat recovered from 50-million-year-old limestone cave sedi- 
ments at Vielase, in the Quercy Phosphorites, southwestern 
France. The fossil fauna from the Vielase deposit is biocorrelated 
as late early Eocene in age, European reference level MP10,*° °° 
which is ~2 million years younger than the Green River bats and 
~1 million years older than those from Messel.°° It is the first 
stem bat described from cave (rather than lake, river, or estua- 
rine) sediments, a depositional environment that has resulted in 
three-dimensional (3D), unflattened preservation of its remains. 
The new bat is described on the basis of a damaged but 
uncrushed skull preserving the petrosal, as well as a rostrum, 
premaxilla, maxillae, dentaries, isolated teeth, and skeletal ele- 
ments. Morphological and total evidence phylogenetic analyses 
consistently place the new taxon outside modern (= crown 
group) bats, with our dated total evidence analysis (the result 
of which we consider to be the most robust°“"°*) placing it sister 
to the crown bats. This fossil taxon provides significant new in- 
formation about the echolocation capabilities of an early Eocene 
stem bat and the early evolutionary history of modern bat 
echolocation. 


RESULTS AND DISCUSSION 


Systematic paleontology 
Order Chiroptera Blumenbach, 1779 

Family indeterminate 

Vielasia gen. nov. 

Type species. Vielasia sigei sp. nov. 

Generic diagnosis. As for the type species until others are 
recognized. 

Generic etymology. From the type locality Vielase. 

Geological setting. Vielase is among the oldest fossil deposits 
identified in the karstic terrane of the Quercy Phosphorites, 
southwestern France.°°'°°°” The Vielase fossils were extracted 
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from a russet-colored bone breccia and a light-colored limestone 
by acid processing.”™®® More than 400 bat specimens referable 
to Vielasia sigei were recovered, with the most complete speci- 
mens (including the holotype cranium UM-VIE-250) being ex- 
tracted from the light-colored limestone (B. Marandat, personal 
communication). Taphonomically, these bat fossils are consis- 
tent with in situ accumulation in a cave;°?°° the remains are 
well preserved but dissociated and include 3D-preserved cranial 
and postcranial specimens representing multiple individuals, 
with teeth and bones unworn and without rounding or size filtering 
indicative of transport, nor evidence of digestion or accumulation 
by (e.g.) predatory birds. Much rarer rodent, primate, and carniv- 
orous mammal fossils recovered from the site°° have been used 
to biocorrelate the deposit as reference level MP10 in the Euro- 
pean Palaeogene mammal biostratigraphy.*? 5? 

Vielasia sigei sp. nov. [urn:lsid:zoobank.org:pub:C9BOAACD- 
3581 -45A3-BB31-36505D03C788] 

Synonymy. Archaeonycteridae gen. et sp. indet. in Sigé°° and 
in Legendre et al. (see figure 7)°° 

Type locality and age. Vielase locality, Quercy Phosphorites, 
southwestern France; late Ypresian, early Eocene (MP10, ~50- 
49 mya).°°°°S 

Holotype. Université de Montpellier (formerly Université de 
Montpellier 2), UM-VIE-250, near-complete cranium with left 
P3-MB8, right P4—M3, root of LC1, and right petrosal; length 
18.29 mm (Figures 1, 2, 3, 4, and S1; Table 1; Data S1A and S1F). 

Paratypes and referred specimens (Figures 1, 2, 4, and S1; Ta- 
ble 1; Data S1A). Paratypes: UM-VIE-251, rostrum with right (R) 
C1-M3 and left (L) P3-M2; UM-VIE-252, right maxilla with C1, 
P3-M3, and alveolus for P2; UM-VIE-668, left premaxilla with 11- 
2; UM-VIE-254, partial left dentary with c1, p3-m3, and alveoli 
for i1-3 and p2; UM-VIE-255, left dentary with p3-m3 and 
ascending ramus; UM-VIE-614, left p2; UM-VIE-669, left petrosal; 
UM-VIE-670, left petrosal; UM-VIE-671, left scapula fragment; 
UM-VIE-672, right proximal humerus; UM-VIE-673, left distal hu- 
merus; UM-VIE-674, left proximal radius; VIE-675, left proximal fe- 
mur; UM-VIE-676, right distal radius; UM-VIE-677, left proximal fe- 
mur. Additionally, UM-VIE-17, right dP4; UM-VIE-02, right dp4; 
UM-VIE-20, right m1 or m2; UM-VIE-19, left M1 (figured by 
Sigé,°° see figures 1-4 therein). Referred specimens: UM-VIE- 
270-UM-VIE-566, UM-VIE-605, UM-VIE-608, and isolated teeth 
(see STAR Methods and Data S1A). Paratype and referred speci- 
mens are known only from the type locality. 

Specific etymology. In honor of French paleontologist Dr Ber- 
nard Sigé, in recognition of his outstanding contribution to the 
study of fossil bats, including the Vielase bat described here. 

Diagnosis. Cranium small (length without prexmaxilla = 
18.29 mm), dorsoventrally flat profile, and without anteriorly in- 
flated braincase; premaxilla proclivous, loosely attached (un- 
fused) to maxilla, nasal process present, and palatine process 
not well developed; weak but present anterolateral rostral infla- 
tions, supraorbital ridges, and postorbital processes; supraor- 
bital foramina opening into two anteriorly directed grooves; 
sagittal crest low anteriorly and reaching maximum height at 
posteriormost point; lambdoidal crest distinct; rostrum narrows 
anteriorly from level of C1-P4; large infraorbital foramen dorsal 
to P4; orbit large and deeply excavated anteriorly; basisphenoid 
pit(s) absent; petrosal loosely attached, cochlea relatively 
expanded, ratio of cochlear width to basicranial width ~0.25, 
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Figure 1. Vielasia sigei n. gen. et sp. 

A-D) UM-VIE-250, holotype, ventral, dorsal, lateral, and petrosal region. 
E) UM-VIE-251, rostrum. 

F) UM-VIE-668, premaxilla with upper incisors. 

G and H) UM-VIE-254, dentary, buccal, and occlusal views. 

I) UM-VIE-255, dentary and lingual view. 

J) Indicative skeletal elements represented in the fossil sample (orange); letters refer to figured specimens. 
K) UM-VIE-671, scapula, glenoid fossa. 

L and M) UM-VIE-672, proximal humerus. 

N) UM-VIE-673, distal humerus. 

O) UM-VIE-674, proximal radius. 

P) UM-VIE-275, proximal femur. 

All scale bars, 5 mm. 

Diagram by N.P. Archer. 

See also and 
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Figure 2. Phylogenetic relationships of Vielasia sigei among bats 


Consensus tree of post-burn-in trees from Bayesian-dated total evidence analysis of 2,665 craniodental characters plus 36,860 base pairs of DNA sequence data 


from 27 nuclear loci from O'Leary et al.°' 


Values at nodes represent Bayesian posterior probabilities. 3D model of Vielasia sigei cranium (compiled from UM-VIE- 


250, 251, and 268) and mandible (from UM-VIE-254 and 255), see STAR Methods. 
Abbreviations: EECO, early Eocene climate optimum (53.26-49.14 mya‘); LE, laryngeally echolocating; mya, million years ago. 


Background photo, Causses du Quercy, P. Lasvenes. 
Scale bar, 5 mm. 
See also Figure S2 and Data S1 and S2. 


2.25 turns in the cochlea, cochlear canal decipherable through 
petrosal bone, and secondary bony lamina extending almost to 
canal apex; mastoid not completely covered by paroccipital 
and supraoccipital; mandible with tall ascending ramus, anterior 
margin of coronoid process nearly vertical, and height at least 
twice tooth crown height; angular process at level of tooth al- 
veoli, digit-shaped, and slightly flared laterally; condylar process 
well dorsal to tooth row; horizontal ramus uniform in depth; ante- 
rior mental foramen ventral to incisors, very large, and round; 
lateral mental foramen large and ventral to p2 or c1-p2; and 
mandibular foramen ventral to tooth row, elongate. 

Relatively complete, unspecialized dentition, with molariform 
dP4/4. Adult dental formula 12/3 C1/1 P2/3 M3/3; 11-2 tall, sti- 
letto-like, without cingula, and weakly bilobed; C1 anteropos- 
teriorly elongate, with single vertical cusp, rounded buccal 
and subplanar lingual faces separated by crests, cingulum dis- 
continuous, and attenuated buccally; P2 usually absent (peg- 
like if present, as in UM-VIE-252); P3 large, with single cusp, 
and three roots; P4 with three roots, wider than long, 
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T-shaped in occlusal view, with tall paracone, and low 
expanded lingual lobe with small protocone; M1-M2 wider 
than long, with large protocone, lingually expanded heel, 
deep and wide single ectoflexus between parastyle and meta- 
style, lacking cuspidate mesostyle and hypocone, and paral- 
oph and metaloph variably present; M3 with paracone, small 
metacone and protocone, three buccal crests only 
(lacking postmetacrista), and similar in width to M2; lower inci- 
sors small, tri-lobed; c1 with posteriorly extended heel, 
cingulum well-developed lingually, absent buccally, and small 
cingular cuspid posteriorly; p4 > p3 > p2; p4 with small but 
distinct metaconid and paraconid, with tall entocristid-like crest 
from metaconid to posterolingual cingulid and cuspid, and 
short talonid with narrow groove or basin perpendicular to 
tooth row; m1-m3 nyctalodont (postcristid links hypoconid 
directly to hypoconulid, isolating entoconid), talonid cusps all 
cuspidate, trigonid cusps taller than talonid cusps, protoconid 
and metaconid similar in height, and trigonid and talonid similar 
in width; cristid obliqua joins postvallid (posterior face of 
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Figure 3. Inner ear attributes in Vielasia sigei compared with extant bats and other mammals 

(A) Virtual models of Vielasia sigei (clockwise from top left): ventral cranium with right petrosal in situ; basilar membrane path within cochlea; endocast of right bony 
labyrinth; and CT slice through modiolar section of cochlea showing foramina (arrows) in bony wall of Rosenthal’s canal. 

(B) Log basilar membrane length (mm) vs. log body mass (g) in extant mammals, with symbol size proportional to the degree of cochlear coiling; phylogenetic 
generalized least squares (PGLS) regression line with 95% prediction intervals for mammals (Data S1G; see also Figure S3A). 

Color key for (B)-(D): aqua, non-laryngeally emitting bats (pteropodids); dark blue, laryngeally emitting bats (triangles, nasal emitters; inverted triangles, oral 
emitters); tan, high-frequency-hearing cetaceans; yellow, low-frequency-hearing cetaceans; black, other mammals; pink diamond, Vielasia sigei. 

(C) Semicircular canal radius of curvature and cochlea size in extant bats (data from Davies et al.°°) and in Vielasia (STAR Methods; Data S1M). Box plots. ASC, 
anterior semicircular canal; LSC, lateral semicircular canal; PSC, posterior semicircular canal; cochlea size.°* 

(D) Plot of mammal species including bats on first two axes of linear discriminant analysis (LDA) of PGLS regression residuals of each semicircular canal log size 
against cochlear log size; percentages represent average between-groups variance explained (see also Figure S3B). 

Mammal silhouettes from PhyloPic (http://www.phylopic.org): pteropodid https:/Awww.phylopic.org/images/1dda6a3d-e407-47fb-add8-fb6b0572f688/pteropus- 
medius by Margot Michaud, CCO 1.0 Universal Public Domain Dedication; yangochiropteran https://www.phylopic.org/images/1 8bfd2fc-f184-4c3a-b51 1- 
796aafcc70f6/corynorhinus-townsendii by Yan Wong, CCO 1.0 Universal Public Domain Dedication; mouse http://www.phylopic.org/image/92989e35- 
4e68-4a2d-b3a2-191ba9da671a/, Mus musculus, by Kamil S. Jaron, Public Domain Dedication 1.0 license; elephant http://www.phylopic.org/image/ 
7c9ab1 82-1 75d-4f02-96d0-09c1e5212bfi/, Elephas maximus, by T. Michael Keesey, Creative Commons Attribution 3.0 Unported license; whale (http://www. 
phylopic.org/image/374accb5-1 6d5-4cb9-a6/a-e881ddfeci 14/, Balaenoptera acutorostrata, by T. Michael Keesey, Creative Commons Attribution 3.0 Unported 
license; dolphin http://www.phylopic.org/image/Scaf4fbod-ca3a-48b4-925a-50fbe9acd887/, Delphinus capensis, by T. Michael Keesey, Creative Commons Attri- 
bution 3.0 Unported license. 

See also Data S1. 
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Figure 4. Skull and mandible shape in Vielasia sigei compared with extant bats 

(A) Cranial morphospaces of bats based on phylogenetic principal component analysis (pPCA) score (phylogenetic principal component [pPC] 1 vs. pPC2) and 
position of Vielasia (red diamond) in morphospace with respect to echolocation call emission type: non-LE, no laryngeal emission; LE-nasal, nasal laryngeal- 
emitters; LE-oral, oral laryngeal-emitters; squares, yangochiropterans; triangles, LE yinpterochiropterans; inverted triangles, non-LE yinpterochiropterans. 
Reconstructed composite cranial model of Vielasia sigei, colors indicate the three virtually assembled elements UM-VIE-250, UM-VIE-251, and UM-VIE-668. 
(B) Mandibular morphospace of bats based on pPCA score (pPC1 vs. pPC2) and position of Vielasia (red diamond) in morphospace with respect to diet. Re- 
constructed composite mandible model of Vielasia sigei, colors indicate the two virtually assembled elements UM-VIE-254 and UM-VIE-255 (and mirror). 


Comparative data from Arbour et al.°° 
See also STAR Methods, Figure S5, and Data S1. 


trigonid) at center; m1—m2 with complete buccal cingulid; and 
m3 talonid reduced but hypoconulid distinct. 

Postcranium relatively unspecialized among bats. Glenoid 
fossa of scapula with only rudimentary dorsal articulation facet 
for greater tuberosity of humerus; proximal humerus with hemi- 
spherical head, greater tuberosity extending just proximal to 
head, shallow fossa for supraglenoid tuberosity of scapula, distal 
humerus articular surface only slightly offset with respect to 
shaft, central capitulum spherical, epitrochlea broad, distal 
spinous process separate from trochlea, not reaching distal 
articular surface; proximal radius epiphyseal surface triangular, 
rounded apices, and without conspicuous ridge separating artic- 
ular facets for capitulum and trochlea; head of proximal femur 
large, spherical, directed anteromedially, with fovea capitis and 
neck, trochanters large, widely spread, lesser trochanter not ex- 
tending proximally beyond base of head, greater trochanter ex- 
tending to half head height, proximal shaft straight, and slightly 
flattened. 

For differential diagnosis and expanded description, see STAR 
Methods and Figure S1. 

Body mass estimate. Long bone shaft diameter provides the 
most accurate estimate of body mass in bats.°° °° In the absence 
of exact-midshaft measurements for Vielasia sigei, we used the 
proxies of lower last premolar (p4), first molar (m1), and upper first 
molar (M1) areas°° (18.31, 17.39, and 20.04 g, respectively), as 
well as humerus minimum and maximum diameter closest to 
mid-shaft (1.44 and 1.70 mm, giving mass estimates of 16.09 
and 23.71 g, respectively). The geometric mean of these proxies 
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indicates a body mass of 18.94 g for Vielasia. This compares with 
the median value of 13.8 g for 905 extant bat species and a range 
of 9-90 g for Eocene stem bats.°’ ©° Among stem bats, Vielasia 
sigei is smaller in skull and toothrow length (18.3 and 7.9 mm, 
respectively; Table 1) than Icaronycteris index (21.2 and 
8.3 mm, see tables 1 and 2 of Jepsen®’) but larger than Palaeo- 
chiropteryx tupaiodon (12.5-17 and 7.5 mm, see table 3 of Rus- 
sell and Sigé’°). Body mass estimates based on long bone shaft 
diameter for the two latter species are 24-27 and 9-16 g, respec- 
tively,°’ and our intermediate estimate of 18.94 g for Vielasia sigei 
is not inconsistent with these reconstructions. 


Phylogenetic relationships of Vielasia sigei and identity 
as a stem bat 

To formally test the relationship of Vielasia sigei n. gen. et sp. to 
crown Chiroptera, we used a matrix of 2,665 craniodental charac- 
ters (925 parsimony informative, 433 variable but non-parsimony 
informative; MorphoBank Project 4407, http://morphobank.org/ 
permalink/?P4407), based on that of O’Leary et al.,°' that includes 
V. sigei and 21 other bat ingroup terminals (6 extant taxa represent- 
ing the major extant lineages within the crown clade, plus 15 
Eocene fossil taxa), plus 8 non-bat eutherian outgroup terminals 
(3 extant, 5 fossil; see STAR Methods). Vielasia sigei could be 
meaningfully scored (i.e., excluding unknown and inapplicable 
characters) for 1,421 characters, rendering it 53.3% complete 
and hence the most complete of all the fossil bats included in our 
matrix (Data S1B), with the others ranging in completeness from 
45.9% (Icaronycteris index) to 5.7% (Microchiropteryx foliae). 


Please cite this article in press as: Hand et al., A 50-million-year-old, three-dimensionally preserved bat skull supports an early origin for modern echo- 


location, Current Biology (2023), https://doi.org/10.1016/j.cub.2023.09.043 


Current Biology 


© Cellress 


OPEN ACCESS 


Table 1. Measurements (mm) for holotype and paratypes of Vielasia sigei n. gen. et sp 


Mandibular 
Cranial parameters UM-VIE-250 UM-VIE-251 UM-VIE-252 Parameters UM-VIE-254 UM-VIE-255 
Cranium L 18.29 - - dentary L - (15.24)a 
Cranium W (mast) 9.01 - - dentary D 2.108 2.131 
Cranium W (zygmax) 12.30 - - dentary H - 5.950 
Cranium W (zygmin) 10.0 - - Rc1-m3 8.816 - 
Interorbital W 4.24 - - Rp4—m3 5.888 5.762 
Rostrum W 7.62 - - ciL 1.172 - 
Rostrum H (M2) 6.34 - - ciW 0.895 - 
Braincase H 6.65 - - p2L - - 
RC1-M3 (7.90) - - p2W - - 
RP4-M3 5.30 - - p3L 1.114 1.097 
RP3L 1.03 1.028 - p3Ww 0.691 0.706 
RP3W 0.798 0.810 - p4L 1.210 1.276 
RP4L 1.56 - - p4w 0.971 0.880 
RP4W 1.640 - - miL 1.680 1.735 
RM1L 1.71 1.730 - m1iWT 1.521 1.525 
RM1W 2.342 2.387 - miWt 1.316 1.343 
RM2L 1.63 1.687 - m2L 1.624 1.676 
RM2W 2.612 2.571 - m2WT (1.335) 1.523 
RM3L 1.0 - - m2Wt 1.052 1.060 
RM3W 2.284 - - m3L 1.589 1.575 
LC1-M3 - 7.860 (7.937) m3WT 1.065 1.018 
LP4-M3 5.28 5.370 5.448 m3Wt 0.749 0.786 
LC1L - 1.45 1.495 - - - 
LC1W - 0.977 1.033 - - - 
LP3L - 1.018 1.060 - - - 
LP3W - 0.808 0.862 - - - 
LP4L 1.53 1.649 1.635 - - - 
LP4W 1.637 1.558 1.616 - - - 
LM1L 1.71 1.741 1.720 - - - 
LM1W 2.381 2.429 2.359 - - - 
LM2L 1.64 1.676 1.685 - - - 
LM2W 2.622 2.512 2.554 - - - 
LM3L 1.0 1.032 1.040 - - - 
LM3W 2.258 2.294 2.393 - - - 


Abbreviations: C, canine; P, premolar; M, molar; upper case, upper toothrow; lower case, lower toothrow; mast, across mastoid; L, length; W, width; 
WT, trigonid width; wt, talonid width; zyg, width across zygomatic arch; numbers in parentheses indicate estimates. 


“Estimate from aligned meshes of UM-VIE-254 and UM-VIE-255. 


Given the widely recognized beneficial impact in phylogenetics of 
combining and analyzing morphological and molecular data simul- 
taneously,’'’° and the extensive DNA sequence data available 
for our extant terminals, we used a total evidence approach, 
combining our craniodental matrix with 36.9 kb of DNA sequence 
data from 27 nuclear loci, taken directly from O’Leary et al.°' The 
inclusion of temporal information has also been shown to increase 
phylogenetic accuracy,°* and so our focal analysis was a dated 
Bayesian analysis of our total evidence dataset using a fossilized 
birth-death model implemented in MrBayes,’”’® with all tips and 
selected nodes assigned age priors and with separate Indepen- 
dent Gamma Rates clock models applied to the morphological 


and molecular partitions (see STAR Methods). The topology that 
resulted from this analysis is shown in Figure 2. 

In this analysis, Vielasia was relatively weakly placed (Bayesian 
posterior probability [BPP] = 0.63) as sister to crown bats, but 
there was much stronger support (BPP = 0.95)—the strongest 
of any grouping outside the crown clade—for a clade comprising 
crown bats, Vielasia and Cambaya (Figure 2). However, there 
was also maximal support (BPP = 1.00) for monophyly of the 
crown clade, excluding Vielasia and Cambaya. Thus, Vielasia is 
placed as a close relative of, but outside, crown Chiroptera. 
The remaining fossil bats formed a paraphyletic assemblage 
outside this clade, with Onychonycteridae as the first family to 
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diverge, which is congruent with previous suggestions.?:2°°” 


Several relationships outside the crown clade are weakly sup- 
ported (BPP < 0.75), as expected given that they involve fossil 
taxa represented by morphological data only,’°°? and particu- 
larly because most of these fossil taxa are known from dental re- 
mains only and hence suffer from high incompleteness (see Data 
S1B); a similar result was found in another recent study of fossil 
bats.°° Thus, although this analysis provides strong support for 
the position of Vielasia within Chiroptera, confident resolution 
of relationships among other stem bats will likely require the dis- 
covery of more complete specimens. 

Assuming the topology presented in Figure 2, monophyly of 
crown Chiroptera to the exclusion of Vielasia is supported by 
49 unambiguous craniodental synapomorphies (Data S2), and 
topologies in which Vielasia was constrained to fall within 
crown Chiroptera (as either a yinpterochiropteran or a yango- 
chiropteran) were rejected by Bayes factors calculated using 
stepping stone analysis (2In(Bio) > 200; Data S1E). In addition 
to this focal analysis, we also carried out maximum parsimony 
(MP) analysis of the full craniodental dataset and MP analysis of 
a “homoplasy-filtered” version of the craniodental dataset in 
which characters that show homoplasy relative to the current 
consensus view of extant bat phylogeny (i.e., with monophy- 
letic Yangochiroptera and Yinpterochiroptera; e.g., Shi and 
Rabosky®* and Amador et al.°*) were excluded, and undated 
Bayesian analysis of the total evidence dataset (see STAR 
Methods and Figure S2). These subsidiary analyses also sup- 
port Vielasia as a stem bat, and a constrained position of Vie- 
lasia within the crown clade is rejected by Templeton tests 
(for the MP analyses) and Bayes factors (for the undated 
Bayesian analysis; see STAR Methods; Data S1E). 

The age of crown bats was poorly constrained in our dated to- 
tal evidence analysis (Figure 2), stretching from the earliest 
Eocene to the late middle Eocene (95% highest probability den- 
sity [HPD] interval: 53.1-38.5 mya). The divergence between Vie- 
lasia and crown bats was estimated as probably occurring during 
the early Eocene (95% HPD interval: 56.5-49.5 mya), and the 
median branch length of 3.8 mya leading to Vielasia indicates 
that it is not a plausible direct ancestor of crown bats. The deep- 
est divergence within Chiroptera, between Onychonycteridae 
and the other taxa, was estimated as probably occurring within 
the Paleocene (95% HPD interval: 66.9-58.9 mya), although chi- 
ropterans are as yet unknown from this epoch. 


Postcranial morphology and flight capabilities 

Like all other Eocene stem bats so far known, Vielasia sigei ap- 
pears to have been capable of powered flight: its skeletal ele- 
ments exhibit specializations for flight found in all bats, extant 
and extinct,*°*° although its postcranium is relatively unspecial- 
ized among bats, being more similar to species of /caronycteris 
and Archaeonycteris than to extant bats (see STAR Methods). 
In Vielasia sigei, several skeletal features are indicative of consid- 
erable rotational movement in the wing joints, with few special- 
ized locking mechanisms, suggesting that it was capable of rela- 
tively maneuverable flight.°°°’ These features include, in the 
shoulder joint, a hemispherical head on the proximal humerus 
and scapula with correspondingly simple, pear-shaped glenoid 
cavity, and, in the elbow joint, a distal humerus with spherical 
central capitulum and broad medial epicondyle with low spinous 
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process, and proximal radius epiphyseal surface without con- 
spicuous ridge separating articular facets for capitulum and 
trochlea (Figures 1K-1P; STAR Methods). 


Petrosal morphology and echolocation capabilities 
Vielasia’s petrosal is relatively isolated from the rest of the skull 
(Figures 1A and 1D), with the enlarged pars cochlearis contact- 
ing the basicranium via four tiny bony splints and otherwise pre- 
sumably loosely attached to the basicranium via ligaments, as in 
all extant bats known to echolocate laryngeally.?°* In non-echo- 
locating pteropodids and Onychonycteris, the petrosal is 
instead fused or sutured to the basicranium, as in other euthe- 
rians.*'7:° Cochlear width (oblique diameter along second 
half-turn of coiled cochlea duct**; 2.50 mm; Data S/F) relative 
to basicranial width (across mastoids**; 9.01 mm; Data S1F) is 
greater in Vielasia than in non-echolocating pteropodids, and 
falls within the range for extant laryngeally echolocating (LE) 
bats (e.g., vespertilionids, emballonurids, rhinopomatids, some 
phyllostomids, and some hipposiderids) ( see figures s8 and s9 
of Simmons et al.°, and figure 12.5.15 of Habersetzer et al.*’). 
It is greater than in Onychonycteris, Icaronycteris, and Archaeo- 
nycteris spp., and more similar to that reported for Hassianycteris 
and Palaeochiropteryx species (see figure 12.5.15 of Haberset- 
zer et al.*’). Overall cochlea size (calculated as the average of 
the diameter of the first cochlear turn, the second cochlear 
turn, and slant height, following Spoor et al.°* and Davies 
et al.°°; Data S1F and S1M) is within the interquartile range for 
extant LE bats but below the median value for all bats 
(Figure 3C). 

In a virtual endocast of Vielasia’s petrosal (Figures 3A and 
S1F), the coiled cochlea represents 74% of the volume of the 
labyrinth (Data S1F; see below). The number of turns in the co- 
chlea is 2.25 (Data S1F), a degree of coiling that is found in mem- 
bers of both groups of extant LE bats, namely Yangochiroptera 
and Yinpterochiroptera; it is greater than in many non-LE ptero- 
podids and some LE carnivorous megadermatids and nycterids 
(see suppl. table s2 of Davies et al.'°, Habersetzer et al.*”, and 
Ekdale®°). Vielasia lacks the greatly enlarged cochlear basal 
turn found in extant bats that use specialized constant frequency 
(CF) calls (e.g., rhinolophids®°). In Vielasia, the secondary bony 
lamina in the cochlear canal extends well beyond the basal 
turn, and almost to the apex, as it does in extant LE bats but 
not in pteropodids.*' 

Mapping Vielasia and extant LE and non-LE bats onto a plot of 
log basilar membrane length against log body mass for placental 
mammals (Figure 3B; data from Davies et al.'°) shows Vielasia 
lying above the mammalian regression line, as do most LE 
bats. That is, compared with non-echolocating placentals, 
basilar membrane length tends to be longer than expected for 
body mass in most LE bats and in Vielasia. Supporting this, we 
found the intercepts of PGLS (phylogenetic generalized least 
squares) regression lines of bats and of highly specialized mam- 
mals (LE bats + high-frequency- and low-frequency-hearing ce- 
taceans) to slightly, but significantly, differ from those of non- 
specialized mammals, whereas slopes do not differ (Figure S3A; 
STAR Methods; Data S1H). 

Cochlear coiling and basilar membrane length are correlated 
with echolocation call parameters and hearing frequencies. '° Us- 
ing the regression equations of Davies et al. '° for extant bats, and 
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based on cochlea coiling (2.25 turns) and basilar membrane 
length (12.38 mm) in the fossil (Data S1F and S1G), we estimate 
the call frequency range for Vielasia sigei to be ~30-56 kHz 
with peak energy of ~42 kHz. These call frequencies are within 
the range recorded by Davies et al. '° for both extant yangochir- 
opterans and LE yinpterochiropterans (see supplemental table 
s8 of Davies et al.'°) and overlap the range those authors recon- 
structed for the hypothetical ancestor of modern bats (~40-65 
kHz with peak energy ~50 kHz). In terms of hearing, we estimate 
a high-frequency hearing limit for Vielasia sigei of ~60-105 kHz (at 
30 and 60 dB, respectively) and low-frequency hearing limit of 
~0.5-1.5 kHz (at 60 and 30 dB, respectively) (Data S1K-S1L; Fig- 
ure S4), given our body mass estimate of 18.94 g. 

Other neuroanatomical properties related to hearing and 
echolocation that can be assessed in Vielasia’s exceptionally 
well-preserved inner ear include the connection between the 
spiral ganglion (housed in Rosenthal’s canal) and the cochlear 
nerve.°' Vielasia sigei exhibits a “foraminal wall” characterized 
by many small apertures piercing the bony wall of the canal 
(Figure 3A), a morphology regarded to be the plesiomorphic 
condition in bats.°' Vielasia sigei shares this cochlear 
morphology with LE yinpterochiropterans, pteropodids, and 
other mammals, but not yangochiropterans, in which the perfo- 
rations in the bony wall are larger (“fenestral wall”) and/or the 
wall is absent (“wall-less”); in yangochiropterans, combinations 
of these two derived states along the canal are thought to 
contribute to the high diversity of echolocation strategies 
seen in those bats.°! 

In Vielasia the petrosal’s delicate vestibular organ is equally 
well preserved and can be precisely measured, a unique occur- 
rence in an early Eocene bat fossil. As occurs in all extant LE bats 
but not in non-LE bats,°° the three semicircular canals deviate 
from circularity in Vielasia (Figures 3A and S1F; STAR Methods; 
Data S1F). Davies et al.°° also found that in LE bats the cochlea 
was relatively much larger than the semicircular canals, as 
measured by the radius of their curvature (R).°* Among mam- 
mals, Vielasia sigei groups among LE bats according to PGLS 
regression residuals of semicircular canal size vs. cochlea size 
(Figure 3D), and when phylogenetic signal is accounted for in 
the analysis, little separation of LE and non-LE mammalian 
groups occurs (Figure S3B; Data S11 and S1J). The anterior 
semicircular canal is larger than that of the lateral and posterior 
canals (Figure 3C), as it is in all yangochiropterans and pteropo- 
dids (Data S1M), rather than the lateral canal being typically 
largest (as in rhinolophids/CF echolocators).°7 

Eye size in Vielasia, estimated here to be less than or equal to 
~2.3 mm diameter (reconstructed using the method of Brooke 
et al.°°; STAR Methods), is smaller than the minimum eye size 
evidently required for aerial pursuit of night-flying insects by extant 
bats without the assistance of echolocation.”° It is also smaller 
than the eye size predicted for the common ancestor of bats 
(3.13 mm).” Relative to body size, Vielasia’s estimated eye size 
is smaller than the smallest modern pteropodid eye, relatively 
larger than in bats that use CF calls (e.g., rhinolophids and hippo- 
siderids) or DH calls (frequency-modulated, dominant harmonic 
calls, e.g., vespertilionids and molossids), and similar to many 
extant bats that use multiharmonic (MH) calls (e.g., emballonurids 
and phyllostomids) (see figure 5, supplement 1 of Thiagavel 
et al.°). 
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Echolocation emission and diet 
3D geometric morphometric (GM) analysis of Vielasia sigei’s cra- 
nial shape, compared with 202 extant bat species of known echo- 
location type® (Data S1N and S1P), suggests V. sigei emitted 
its echolocation call through its mouth rather than its nose 
(Figures 4A and S5A) based on the two first phylogenetic principal 
component (pPC) axes, which largely reflect cranial elongation 
(pPC1, 37.1%) and rostral dorsoventral flexion (pPC2, 23.1%). 
Of the 202 extant bat species compared, Vielasia lies closest in 
pPC1-pPC2 morphospace to the vespertilionid Myotis keenii, 
the molossid Mops condylurus, and the natalid Natalus strami- 
neus, all of which are oral emitters.°°'°*°° As in extant oral-emit- 
ting bats (e.g., vespertilionids, emballonurids, and mormoopids), 
the hard palate in Vielasia is approximately aligned with the basi- 
cranium (small positive tilt of the rostrum on the basicranial axis), 
the foramen magnum is directed posteriorly, and the plane of the 
lateral semicircular canal is relatively parallel with that of the ba- 
sicranium.°“"°° Via this arrangement, orofacial and cranial struc- 
tures are aligned, such that undistorted echolocation calls are 
emitted from the open mouth directly in the line of flight.°° Vielasia 
lacks the ventral flexion of the rostrum characterizing Old World 
nasal emitters (e.g., rhinolophids and hipposiderids) and New 
World nasal emitters (e.g., phyllostomids).°°24"°° 

Arbour et al.°° found that mandible shape better reflects diet in 
extant bats than cranial shape. Compared with 191 extant bat 
taxa (Data S10 and S1P), Vielasia groups in mandibular shape 
morphospace primarily with insectivorous bat species (Figures 4B 
and S5B). Vielasia is intermediate in anteroposterior elongation of 
the dentary (pPC1, 35.3%) and mandible body and coronoid 
height (pPC2, 16.8%), grouping with insectivorous, carnivorous, 
and frugivorous bats (Figure 4B). However, its intermediate dorso- 
ventral mandible flexure (pPC3, 11.9%) and tooth row length 
(pPC4, 11.6%) place it within the morphospace occupied by insec- 
tivorous bats only (Figure S5B). Stem bats whose gut contents are 
preserved in other fossil Lagerstatten also appear to have been pri- 
marily insectivorous, with Messel bats such as Palaeochiropteryx 
spp. evidently consuming predominantly small moths and caddis- 
flies*°*” and the larger Archaeonycteris trigonodon and Hassia- 
nycteris spp. interpreted to have included larger moths and beetles 
in their diet.4°4°* 


Cave dwelling in a changing world 
Bats first appeared in the fossil record directly following the 
Paleocene/Eocene Thermal Maximum (PETM) at 55.93 mya.°° 
They radiated during the Early Eocene Climatic Optimum 
(EECO, 53.26-49.14 mya®?; Figure 2) during global “hothouse” 
conditions,°° and rapidly appeared on most continents.”° These 
early Eocene bats are typically fossilized in fluvio-lacustrine or 
estuarine deposits and are interpreted to have roosted in 
trees,27°146°7°5 in the paratropical Eocene forests that once 
stretched to the poles.°*°° 

Vielasia sigei, however, demonstrates that some bats had 
moved into caves by the late early Eocene. Vielasia is the oldest 
bat recovered from a cave deposit, with cave dwelling otherwise 
first recorded in early members of the modern bat family Rhino- 
lophidae in the middle Eocene of China.“ At least 23 individuals 
of V. sigei occur in the Vielase deposit, including at least one ju- 
venile represented by its molariform deciduous posterior denti- 
tion (dP4/4).°° This dental condition is plesiomorphic for 
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placental mammals and contrasts with the condition in all extant 
bats in which the entire deciduous dentition (dl1/1 to dP4/4) is 
uniquely specialized (typically hook-like'°°); the latter enables a 
newborn bat to cling passively but safely to its mother during 
flight and suspended roosting. In Vielasia, wear facets on 
dP4/4°° indicate that juveniles began to hunt insects before 
weaning, as deduced also from similar dental evidence and 
confirmed by gut contents in pups of the Eocene bat Palaeochir- 
opteryx tupaiodon.'°' 

An important evolutionary driver for bats like Vielasia sigei to 
first congregate in caves may have been the transition to globally 
cooler and less stable conditions, which began around 49 mya 
between the early and middle Eocene and was preceded by 
high-amplitude variations in temperature.°* This marked the 
end of the hothouse EECO, although under the cooler green- 
house conditions of the early middle Eocene, the world was still 
much warmer than now.°° 

Alternatively, caves in these middle latitudes could have pro- 
vided attractive new roosts for evolving bats with advanced 
echolocation, flight, and metabolic capabilities, as well as 
changing behavioral and social needs.°° Among benefits for 
extant bats, cave roosting offers not only a thermally stable envi- 
ronment but also protection from predators, increased mating 
opportunities and success in rearing young, foraging, social 
and behavioral information transfer, and reduced competition 
for tree roosts from birds and other mammals. '°%'°? In higher lat- 
itudes today, it facilitates hibernation, reducing energy costs for 
bats, '°* with cave-dwelling species also generally migrating less 
far than tree-roosting species.'°° Cave roosting is correlated 
with increased longevity in bats. '°°'°” Today, caves and under- 
ground habitats are used by almost half of all bat species. '°° 


Early acquisition of echolocation in bats 
Based on morphological features including relative cochlear 
size, shape of the orbicular process of the malleus, and stylohyal 
morphology,” 1840-4" previous paleontological studies have sug- 
gested that Eocene stem bats were capable of laryngeal echolo- 
cation,”''°*°“” with the possible exception of Onychonycteris 
finneyi,°"*° but see Veselka et al.'° Relative cochlear size and 
an articulation between the tympanic bone in the middle ear 
and a cranially expanded stylohyal bone (part of the hyoid arch 
that supports the larynx) remain definitive and uncontested hall- 
marks of laryngeal echolocation in bats.'®“° The stylohyal and 
middle ear bones are not yet known for V. sigei, but its three- 
dimensionally well-preserved cranium and petrosal (cochlea 
and semicircular canals) has enabled this early bat’s inclusion 
in analyses of large datasets of cranial and inner ear attributes 
related to auditory and echolocation capabilities in extant bats. 
In these analyses, V. sigei grouped with LE rather than non-LE 
bats, and there is little to suggest that Vielasia used a type of 
echolocation different from that used by modern LE bats. This 
finding is congruent with previous studies suggesting that 
some non-crown bats may have had modern or near-modern 
echolocation capabilities*’ and contrasts with some proposals 
based on embryological and molecular analyses of extant bats 
that, if stem bats did use echolocation, it was of a primitive 
kind (such as tongue-clicking). 1° 

At 50 mya, Vielasia’s uncrushed, near-complete skull is the 
oldest yet known, and its description here significantly expands 
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empirical data available for inclusion in chiropteran phylogenetic 
analyses, evolutionary modeling of the development of bat echo- 
location, and in embryological and molecular analyses investi- 
gating high-frequency hearing and the inner ear of bats. 


Conclusions 

Our results indicate that all attributes correlated with echoloca- 
tion that are measurable in the stem bat Vielasia sigei fall within 
the range observed in extant laryngeally echolocating bats, 
and this similarity in reconstructed hearing capabilities suggests 
V. sigei was also capable of modern laryngeal echolocation. 
The analyses presented here collectively suggest either that 
advanced echolocation evolved once in the common ancestor 
of extant bats (and was used by stem bats such as Vielasia sigei) 
but was lost in pteropodids, or that advanced echolocation 
evolved independently several times in bats, at least once in 
non-crown bats and twice in extant bat lineages. 

The early Eocene Vielasia sigei was a relatively small (~18.94 g), 
insect-eating, evidently oral-emitting echolocator, with expanded 
cochlea, reduced eye size, and capable of maneuverable flight. 
Our calculations suggest it may have used a multiharmonic echo- 
location call, with an estimated frequency range of ~30-56 kHz. It 
appears to have been a cave-dweller, living in small colonies that 
included semi-independent young. Given that our preferred phylo- 
genetic hypothesis places Vielasia as the closest known relative of 
crown Chiroptera, at least some of these features may also have 
characterized the last common ancestor of modern bats. 
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EXPERIMENTAL MODEL AND SUBJECT DETAILS 


The main subject of this study, the holotype specimen of Vielasia sigei (UM-VIE-250; Figures 1, 2, 3, 4, and S1), is a well preserved but 
damaged cranium containing the petrosal. This fossil specimen, and others assigned here to V. sigei (Figures 1 and S1), were all 
found at the type locality of Vielase, Quercy Phosphorites, southwestern France,°° and were previously referred to Archaeonycter- 
idae gen. et sp. indet.°°°° 


METHOD DETAILS 


Data collection, CT scanning, model reconstruction, and measurements 

Three dimensional data acquisition from fossils was performed using the „CT facilities (EasyTom 150 CT) at the University of Mont- 
pellier (MRI, ISEM) with 90 kV and 80 pA at a resolution of 8.93 um, and a U-CT (Milabs, Utrecht) at UNSW Sydney with 55 kV and 0.17 
mA, ultrafocused setting at a resolution of 35-50 um. Segmentation, visualization, and linear and volumetric measurements were per- 
formed using Avizo ® 9.3 (Thermo Fisher Scientific-FEl) and MorphoDig© at the University of Montpellier, and MIMICS 21.0 (Mate- 
rialise) at UNSW Sydney. MeshLab 2020.12 and Meshmixer (Autodesk Inc.) were used to align meshes of UM-VIE-250 (cranium), 
UM-VIE-251 (rostrum) and UM-VIE-668 (duplicated and mirrored premaxilla) in order to reconstruct a composite virtual cranial model 
for inclusion in our 3D geometric morphometric analysis of bat cranial shape (see below). Meshes of UM-VIE-254 and UM-VIE-255 
were similarly aligned to reconstruct a composite dentary. Measurements of craniodental material were made on the virtual models 
using 3Matic 13.0 (Materialise) and on the fossils using a Wild 5MA stereomicroscope with Wild MMS235 Digital length measuring set 
(accurate to 0.01 mm). Univariate and multivariate analyses of the dental sample were performed in the statistical program PAlae- 
ontological STatistics (PAST 4.05; Data S1A, and see below). 

One binning of the „CT dataset was applied prior to skull and bony labyrinth segmentation, resulting in a dataset with a 17.86 um 
resolution, using ImageJ v1.48; segmentation was performed using Avizo ® 9.3 tools brush, lasso, and magic wand using varying 
threshold to remove sediment; model reconstruction was performed using Avizo ® 9.3, with an unconstrained smoothing value of 
2.75. Measurements of the bony labyrinth (linear and angles) were made on the virtual model using 3Matic (Materialise), following 
the methods of Spoor and Zonneveld, '** Spoor et al.,°* Ekdale,°° and Davies et al.'°°° (Data S1F, and see below); basilar membrane 
length was measured using 86 approximately equidistantly stationary points (as per Davies et al.'°; Figure 3A) placed in Landmark 3.6. 


Comparative material 

In the list below, an asterisk (*) denotes original material examined in addition to (otherwise) primary literature, supplementary pho- 
tographs and/or scan data from MorphoSource and DigiMorph. Families listed alphabetically. List of institutional abbreviations for 
specimen numbers follows the taxon list. 

AEGYPTONYCTERIDAE: Aegyptonycteris knightae Simmons et al.'°° Fayum, Egypt CGM83740*; ARCHAEONYCTERIDAE: Ar- 
chaeonycteris trigonodon Revilliod'*° Grube Messel near Darmstadt, Germany HLMD 1398- Me 33a, b*, HLMD 16069; A. brailloni 
Russell et al.1?7 Mutigny and Avenay, Marne, France MNHN Louis-410 Mu, MNHN Louis-432 Mu*, MNHN Bn-751 Av*, MNHN 
Bn-757 Av* (cast SMF 77/262); A. pollex Storch and Habersetzer®*’ Grube Messel near Darmstadt, Germany SMNK Me 982°; 
A. relicta Harrison and Hooker'*® Creechbarrow, Dorset, England BMNH M35710/12. 7ARCHAEONYCTERIDAE: ?Archaeonycteris 
praecursor Tabuce et al.” Silveirinha Portugal UNLSNC-447*; 2A. storchi Smith et al.'°° Vastan lignite mine, Gujarat, India GU/RSR/ 
VAS 140; Protonycteris gunnelli Smith et al.'°° Vastan lignite mine, Gujarat, India GU/RSR/VAS 436; Xylonycteris stenodon Hand and 
Sigé°° Prémontré, France SLP 1 PR73*, SLP 27 PR184*. EMBALLONURIDAE:Tachypteron franzeni Storch et al.'°° Grube Messel 
near Darmstadt, Germany IRSNB BE 4-119a+b. HASSIANYCTERIDAE: Hassianycteris revilliodi (Russell and Sigé)’° Grube Messel 
near Darmstadt, Germany HLMD 4294 Me 16* (cast SMF 77/269); H. messelensis Smith and Storch'®' Grube Messel near Darmstadt, 
Germany HLMD Me 7480*; H. magna Smith and Storch'*' Grube Messel near Darmstadt, Germany HLMD Me 7539"; H. kumari Smith 
et al.'°° Vastan lignite mine, Gujarat, India GU/RSR/VAS 56, GU/RSR/VAS 561. 7HASSIANYCTERIDAE: Cambaya complexus Bajpai 
et al.'°° Vastan lignite mine, Gujarat, India ITTR/SB/VLM 435. ICARONYCTERIDAE: Icaronycteris index Jepsen?” Green River Forma- 
tion, Wyoming, USA YPM-PU 18150* and UW 21481"; I. gunnelli Rietbergen et al.°° Green River Formation, Wyoming, USA AMNH:F- 
M:145747A,B, ROM:Palaeobiology-Vertebrate Fossils:52666; /. menui Russell et al." Mutigny and Avenay, Marne, France MNHN 
Louis-333*, MNHN Louis-Mu*, MNHN Louis-331*, MNHN Louis-380*, MNHN Louis-384*, MNHN Louis-353*, MNHN Louis-357*, 
MNHN Louis-117*, MNHN Av7201*, MNHN Av6797*, MNHN Av7200*, MNHN Louis-360 Mu*, MNHN Bn-547 Av*, MNHN Bn-340 
Av*, MNHN Bn-319 Av* (SMF 77/259-261). 7ICARONYCTERIDAE: Icaronycteris sigei Smith et al.'?° Vastan lignite mine, Gujarat, 
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India GU/RSR/VAS 137. MIXOPTERYGIDAE: Carcinipteryx maximinensis Maitre et al.'°° from St-Maximin, Gard Phosphorites, 
France, UM-SMXCm A.1.3*, UM-SMXCm A.1.5*. NECROMANTIDAE: Cryptobune thevenini Sigé'** unknown locality, Quercy, 
France ACQ 140*; Necromantis adichaster Weithofer'®° Quercy, France QU16369*, QU16369"; N. gezei Hand et al. '°° Cuzal, Quercy, 
France CUZ-383*; N. marandati Hand et al.'°° La Bouffie, Quercy, France BFI-NmA.1.1*, BFINmA. 1.5*. ONYCHONYCTERIDAE: 
Onychonycteris finneyi Simmons et al.° Green River Formation, Wyoming, USA, ROM 55351A,B*, AMNH FM 142467*; Ageina tobieni 
Russell et al.'°” Mutigny, Marne, France MNHN Mu 5112*, MNHN Louis-481 Mu* (cast SMF 77/263); Eppsinycteris anglica (Forster 
Cooper)'*’ Abbey Wood, London Basin, England BMNH M13776; Honrovits tsuwape Beard et al.'°° Wind River Formation, Wyom- 
ing, USA CM62640*, CM 62641*; ?Honrovits joeli (Smith and Russell)'®° Evere, Belgium IRSNB M 1567 (cast); Marnenycteris 
michauxi Hand et al.°' Pourcy, France UM POY-11*; Volactrix simmonsae Czaplewski et al. Elderberry Canyon, Nevada, USA, 
USNM PAL 544092, USNM PAL 417351. PALAEOCHIROPTERYGIDAE: Palaeochiropteryx tupaiodon Revilliod'*° Grube Messel 
near Darmstadt, Germany, HLMD Me25* (cast SMF 77/271 and 272), Me 266* (cast SMF 77/273 and 274); Palaeochiropteryx spiegeli 
Revilliod'°° Grube Messel near Darmstadt, Germany HLMD Me32 * (cast SMF 77/270); Palaeochiropteryx sambuceus Czaplewski 
et al.°° Elderberry Canyon, Nevada, USA, USNM PAL 417350, USNM PAL 706598. Cecilionycteris prisca Heller'*° Geiseltal Obere 
Mittelkohle, Germany, GMH 3965*; Matthesia germanica Sigé and Russell'*' Geiseltal Obere Mittelkohle, Germany, GMH 3940"; 
Lapichiropteryx xiei Tong“? Tuqiaogou, Shanxi province, China IVPP V10204; Microchiropteryx folieae Smith et al.'°° Vastan lignite 
mine, Gujarat, India GU/RSR/VAS 459; Anatolianycteris insularis Jones et al.°* Pontide, Anatolia, Turkey, EOQU-UCF-7 and 8. 
PHILISIDAE: Dizzya exsultans Sigé'“° Chambi, Tunisia CB 1-15*. VESPERTILIONIDAE: Premonycteris vesper Hand et al.'“* 
Prémontré, France SLP 29 BS 193*, SLP SLP 29 PE 875*, SLP 29 PE 110*; SLP 29 P E 614*; Sonor handae Czaplewski et al.°° Elder- 
berry Canyon, Nevada, USA, USNM PAL 544091. FAMILY INCERTAE SEDIS: cf. Archaeonycteris sp. (sensu Smith et al.2°); cf. Icar- 
onycteris sp. in Russell et al.'*° Meudon, France MNHN 16065 Me’; cf. Archaeonycteridae (Sigé°°; Legendre et al.°°) Vielase, Quercy, 
France VIE 19%, VIE 20*. Australonycteris clarkae Hand et al” Murgon, Queensland, Australia QM F19147, QM F19149*; Jaegeria 
cambayensis Bajpai et al. °° Vastan lignite mine, Gujarat, India IITR/SB/VLM/585; Chiropteran in Tejedor et al. '*° Laguna Fria, Chubut 
Province, Argentina, LIEB-PV 999; Chiroptera indet. 2 in Marandat'“” Fournes, Minervois, France FNR-02*; Eochiroptera in Ravel 
et al.”° El Kohol, Algeria UT-KD-02, UT-KD-03; Altaynycteris aurora Jones et al.°° Junggar Basin, Xinjiang, China IVPP V 27157, 
IVPP V 27158. 


Institutional abbreviations 

AMNH American Museum of Natural History, New York, USA; BMNH The Natural History Museum, London, UK; CGM Cairo Geolog- 
ical Museum, Egypt; CB Chambi collection, Tunisia; CM Carnegie Museum of Natural History, Pittsburgh, USA; EOU, Eskisehir Os- 
mangazi University, Eskişehir, Turkey; FNR Université de Montpellier Il, Fournes Locality Collections, France; GMH Geiseltal 
Museum, Halle, Germany; GPIMUH Geologisch-Palaontologisches Institut und Museum, Universitat Hamburg, Germany; GU/ 
RSR/VAS Garhwal University, Srinigar, India; HLMD Hessisches Landesmuseum, Darmstadt, Germany; IITR/SB/VLM Indian Institute 
of Technology, Roorkee, India; IRSNB M Institut Royal des Sciences Naturelles de Belgique, Brussels, Belgium; IVPP V Institute of 
Vertebrate Paleontology and Paleoanthropology, Beijing, China; LIEB-PV Laboratorio de Investigaciones en Evolución y Bio- 
diversidad, Universidad Nacional de la Patagonia San Juan Bosco, Argentina; MNHN Muséum National d’Histoire Naturelle, Paris, 
France; QM Queensland Museum, Brisbane, Australia; ROM Royal Ontario Museum, Toronto, Canada; SLP, Société Laonnoise de 
Paléontologie collections, Muséum National d’Histoire Naturelle, Paris, France; SMF Me Senckenberg Museum, Messel Department, 
Frankfurt, Germany; SMNK Staatliches Museum fur Naturkunde, Karlsruhe, Germany; UM Université de Montpellier Il, Montpellier, 
France; UNLSNC Universidade Nova de Lisboa, Portugal, Silveirinha new collection; USNM PAL, Department of Paleobiology, Na- 
tional Museum of Natural History; UT University of Tlemcen, Tlemcen, Algeria; UW University of Wyoming Museum, Laramie, Wyom- 
ing, USA; YPM-PU Princeton University (collections now housed in the Yale Peabody Museum, New Haven, CT, USA). 


Differential diagnosis 
Vielasia sigei differs from all other bats in the combination of features noted in the Diagnosis (main text). Morphological features dis- 
tinguishing V. sigei from specific bat taxa are given below. See Figure S1 for cranial, mandibular and dental terminology. 

The dental and postcranial morphology of Vielasia sigei most closely resembles that of early to middle Eocene bats referred to 
several extinct bat families (or currently unassigned to family). Vielasia sigei is distinguished from these groups as follows. 

Differs from onychonycterids (species of Onychonycteris, Honrovits, Ageina, Eppsinycteris, Marnenycteris and Volactrix) in its up- 
per molars being transversely developed (wider than long rather than nearly square); P4 relatively larger; M3 lacking postmetacrista; 
p3 relatively larger, and p2 smaller such that p2<p3<p4; m1-3 nyctalodont, cristid obliqua meeting the postvallid at centre; dentary 
with more vertical anterior margin of coronoid process. Differs additionally from Volactrix simmonsae in p3 having two roots. Differs 
additionally from E. anglica in p3 having only two roots and m1 not smaller than m2. Differs additionally from O. finneyi in 11 and 12 
being of similar size, P2 tiny or absent, P3 relatively larger and with three roots, M1 with ectoflexus, M3 relatively shorter, lower in- 
cisors tricuspid, more expanded cochlea, shoulder joint with rudimentary secondary articulation, and capitulum of distal humerus 
near-aligned with longitudinal axis of shaft (rather than laterally displaced with condyle lateral to lateral ledge of shaft). 

Differs from icaronycterids (/caronycteris spp.) in P2 absent or tiny, P3 relatively smaller, P4 relatively larger, M1-3 lacking hypo- 
cone, and M3 relatively shorter. Differs additionally from /. index in its premaxilla with reduced or absent palatine process, 11-2 not 
recurved lingually, P4 lacking hypocone, dentary coronoid process with more vertical anterior margin and deeper horizontal 
ramus with slightly convex ventral margin, more expanded cochlea, shoulder joint with rudimentary secondary articulation and 
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capitulum/articular surface of distal humerus less laterally displaced. Differs additionally from /. sige in its smaller p4 metaconid. Dif- 
fers additionally from /. gunnelli in lacking post-C1 diastema, and in P3 height less than P4 height. 

Differs from archaeonycterids (species of Archaeonycteris, Protonycteris and Xylonycteris) in P2 absent or tiny, P3 relatively larger 
and with three roots, P4 with more lingually expanded heel, M1-3 with weaker cingulae including around heel, m1-3 with more cuspi- 
date entoconid and hypoconulid, cristid obliqua meeting the postvallid at centre. Differs additionally from Archaeonycteris spp. in its 
more expanded cochlea. Differs additionally from P. gunnelli and X. stenodon in m1-3 nyctalodont, and from X. stenodon in its rela- 
tively wider lower molars. 

Differs from palaeochiropterygids (species of Palaeochiropteryx, Cecilonycteris, Matthesia, Lapichiropteryx, Microchiropteryx, 
Anatolianycteris, and possibly Stehlinia, the latter sometimes referred to Vespertilionoidea”’) in its P4-M3 more transversely devel- 
oped, particularly M2 protocone lobe, 11-2 bilobed, lacking cingulum and cingular cuspule, P2 absent or tiny, P3 relatively smaller, P4 
T-shaped rather than triangular, M1-2 without hypocone development and without double ectoflexi, m1-3 with relatively narrower 
talonids, entocristid and postcristid meet at angle of less than 90°. Differs additionally from L. xiei (and Stehlinia spp.”°) in its 
M1-2 with ectoflexus and metastylar shelf, and from C. prisca in p4 with metaconid. Differs additionally from Stehinia minor in its distal 
humerus with spherical capitulum further laterally offset and spinous process short, not extending distally beyond the articular sur- 
face, and in its broader rostrum, anteriorly convergent [and relatively shorter] upper tooth rows, premaxilla (with palatine process and) 
without medial cleft, anterior palate margin V-shaped, infraorbital formen dorsal to P4 rather than P3 (see figure 12.5.12 of Haber- 
setzer et al.“’). Differs additionally from Anatolianycteris insularis (see figure 2 of Jones et al.°°) in its less elongated p4. 

Differs from hassianycterids (species of Hassianycteris) in its M1-2 with ectoflexus, 11-2 not diamond shaped, P3 with three roots, 
M1-2 with heels less posterolingually developed, M3 with three roots, p4 with metaconid and longer talonid, m1-3 with relatively nar- 
rower talonid and lower entoconid, m3 talonid less reduced, angular process of dentary slimmer, cranium with uninflated braincase. 
Differs additionally from H. messelensis and H. magna in P4, p3 and p4 without exceptionally tall primary cusp, from H. revilliodi in c1 
not recurved, and from H. magna in trigonid and talonid cusps not subequal in height. 

Differs from Cambaya complexus (tentatively referred to Hassianycteridae by Smith et al. '°°) in its shallower dentary without mark- 
edly anteroventrally expanded symphysis, p4 with wider (and longer?) talonid, m1-3 with relatively narrower talonid, m1 with talonid 
narrower than trigonid, angular process of dentary slimmer, inflection point for ascending ramus of dentary ventral (rather than pos- 
terior) to m3, with posteroventral margin rounded. 

Differs from mixopterygids (species of Mixopteryx and Carcinipteryx) in its C1 not recurved, P4 T-shaped, M1-2 more transversely 
developed and without posteriorly directed heel (not quadrate occlusal outline), without protruding mesostyle, postprotocrista 
meeting base of metacone, three lower incisors, three lower premolars, m1-2 without crestiform hypoconulid, m2 talonid relatively 
narrower. 

Differs from aegyptonycterids (Aegyptonycteris knightae) in its M2 being more transversely developed, with deeper ectoflexus, 
lacking ectostyle between parastyle and centrocrista, without large, bulbous hypocone, lacking large metaconule, postprotocrista 
not longitudinal to tooth row, lingual and buccal cingulae present; M3 more reduced and without hypocone shelf. 

Differs from necromantids (species of Necromantis and Cryptobune) in its P3 large and with three roots, P4 wider and shorter, M1-2 
without enlarged heels, centrocrista reaching buccal border, M3 with premetacrista (and metacone), p3 large, p4 short and with 
metaconid, and m1-3 nyctalodont. Differs additionally from N. adichaster in premaxilla and upper incisors present, P2 absent or 
tiny, lower incisors present, gracile cranium, angle between basicranium and facial region not pronounced, supraorbital ridges 
and postorbital processes weak, facial region not with clear pentagonal shield, palate relatively longer with posterior margin posterior 
to toothrow, unexpanded zygoma, paroccipital not overlying mastoid, horizontal ramus of dentary not exceptionally deep, more 
expanded cochlea, and distal humerus with relatively smaller, less laterally displaced, capitulum and longer distal spinous process. 

Differs from the tanzanycterid Tanzanycteris mannardi in its less expanded cochlea (dentition unknown in that species). 

Differs from philisids (Species of Philisis, Dizzya and Witwatia) in its P4-M3 more transversely developed, M1-2 ectoloph without 
subparallel cristae and double mesostyle, p4 with metaconid, and m1-3 nyctalodonrt. 

Differs from bat species unassigned to family as follows. 

Differs from Jaegeria cambayensis (see figure 2 of Smith et al.'7°) in its deeper dentary, p4 with smaller metaconid, m1 talonid nar- 
rower than trigonid. Differs from Chambinycteris pusilli°° in its C1 and c1 relatively taller, P3 large with three roots, M1 without pos- 
teriorly directed heel, p3 large with two roots, p4 with two roots unfused, dentary with deeper horizontal ramus. Differs from Drako- 
nycteris glibzegdouensis”° in its P3 not massive, M1-2 ectoloph without subparallel cristae, M1 with single mesostyle and strong 
parastyle, M2 buccal margin directed anteroposteriorly (rather than obliquely), m1-3 with more cuspidate hypoconulid. Differs 
from Australonycteris clarkae“ in its P4 with parastyle, lingually extended heel/lobe and large protocone/cingular cusp, m1-3 nyc- 
talodont. Differs from Altaynycteris aurora (see figure 2 of Jones et al.°°) in its M1-2 with deep ectoflexus and presence of heel. 

Differs in its m1-3 nyctalodont from early Eocene indeterminate bat taxa such as:?Archaeonycteris sp., Meudon, Paris Basin'*°; 
Eochiroptera indet. 2 (see pl. 5 of Marandat'*’), Fournes, southern France; chiropteran indet. from Laguna Fria, Argentina'“°; and 
Eochiroptera indet. from El Kohol, Algeria (see figures 3-4 of Ravel et al.”°). Differs additionally from the El Kohol eochiropteran its 
M1-3 with smaller paracone and narrower cingulum. 

Vielasia sigei differs from modern bat species in its (archaic) dental morphology. This suite of features includes: molariform dP4 and 
dp4, P4 paracone tall, twice height of M1 paracone, P4 protocone medial to paracone, M1-3 transversely developed (wider than 
long), anteroposteriorly narrow, tall cusps and deep valleys, large protocone and lingually extended heel without hypocone, M1-2 
with buccal and anterolingual cingula, M1 stylar shelf narrow (less than 50% width of tooth), metacone more lingual than paracone, 
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protocone tall (more than 80% height of metacone), postprotocrista extending buccally to metacone base, M2 paracone and meta- 
cone subequal in height, M2 with tall postprotocrista, M3 large, two-thirds area of M1, p3 tall, large (80-40% area of m1), with para- 
conid and tall protoconid, subequal to p4 protoconid height, p4 with paracristid, talonid more than one third trigonid length, with hy- 
poconulid, hypoflexid, small entoconid, m1-3 trigonids becoming more mesiodistally compressed from m1 to m3, m1 with low 
paracristid, only 50-75% height of paracone, metaconid height subequal to protoconid height, talonid subequal in width to trigonid, 
with postentocristid, m2 with with metaconid height subequal to protoconid height, complete but low protocristid, low weak crest 
between hypoconid and hypoconulid, hypoconulid above cingulid level, m3 hypoconulid at or near midline of tooth. 

Vielasia sigei differs additionally from individual modern bat families as follows. 

Differs additionally from pteropodids in its tribosphenic teeth not highly modified (tribosphenic morphology clearly identifiable), 
premaxilla with at least small flange on palatine process (Figure S1G), cranium with expanded cochlea, shoulder joint with rudimen- 
tary second articulation and distal humerus articulation less laterally offset. 

Differs additionally from rhinolophids, hipposiderids and rhinonycterids in its premaxilla with nasal process and two incisors, nasal 
not inflated, cranium with short postorbital process, short p4 with metaconid, dentary with tall ascending ramus, and distal humerus 
with articular surface less laterally offset. Differs additionally from rhinonycterids and hipposiderids in having three lower premolars, 
although Eocene hipposiderid species (e.g., Palaeophyliophora and Pseudorhinolopus) may retain p3. Differs additionally from rhi- 
nolophids in its longer hard palate without deep U shaped anterior emargination, large, two-rooted p3, tall coronoid process of den- 
tary, less expanded cochlea, distal humerus without deep ridges and long spinous process, proximal femur with large trochanters 
and shaft straight. 

Differs additionally from rhinopomatids in its rostrum without anterolateral inflations, cranium with short postorbital process, its two 
upper incisors, C1 relatively larger, two upper premolars, M1-2 without large heel, M3 with metacone, three lower incisors, three 
lower premolars, anterior lower premolar larger, m1-2 talonid narrower, proximal humerus with hemispherical head, second (though 
rudimentary) articulation. 

Differs additionally from craseonycterids in premaxillae not forming complete ring around narial opening, rostrum without antero- 
lateral inflations, braincase not inflated, with postorbital process, less expanded cochlea, two upper incisors, more than one upper 
premolar, M1-2 without large heel, M3 with metacone, three lower incisors, three lower premolars, anterior lower premolar larger, 
m1-2 talonid narrower, dentary with tall ascending ramus, distal humerus articulation less laterally offset and spinous process short. 

Differs additionally from megadermatids in premaxilla with well-developed nasal process and bearing incisors, braincase less 
expanded, more expanded cochlea, C1 without large posterior accessory cusp, M1-2 without large heels, M3 with metacone and 
premetacrista, three lower incisors, three lower premolars, p4 with metaconid, m1-3 without relatively long, open trigonids and small 
talonids, proximal humerus with hemispherical head, distal humerus articulation less laterally offset, proximal femur with well-devel- 
oped trochanters and straight shaft. 

Differs additionally from emballonurids in its cranium with less developed postorbital process, without basisphenoid pits, P3 not 
highly reduced, M1-2 without posteriorly expanded heels, c1 without long posterior heel, three lower premolars, m1-3 talonids rela- 
tively narrow, proximal humerus with hemispherical head, and proximal femur with well-developed trochanters and straight shaft. 
Differs additionally from Taphozous spp. inits two upper incisors, three lower premolars, anterior dentary not abruptly shallow. Differs 
additionally from Tachypteron franzeni in its M3 without heel, three lower incisors, relatively taller c1, larger p2, p4 with metaconid and 
entocristid-like crest on talonid, ventral margin of dentary without marked concavity below incisor region, glenoid fossa of the scapula 
with secondary articular surface, head of proximal humerus hemispherical, greater tuberosity extending proximally beyond head, 
distal humerus with spinous process not reaching level of trochlea. Differs additionally from Pseudovespertiliavus and Vespertiliavus 
spp. in its m1-3 trigonid not anteroposteriorly compressed, and additionally from Vespertiliavus spp. in its three lower premolars, with 
p4 longer than p2. 

Differs additionally from nycterids in premaxilla with nasal process, cranium without basisphenoid depression/pits, more than one 
upper premolar, M1-2 without posteriorly expanded heels, M3 not strongly reduced, cranium lacking deep frontal depression and 
very broad supraorbital ridge extending over orbit, c1 tall, three lower premolars, p4 with metaconid, m1-3 with entoconids not mark- 
edly reduced, m3 length and width less reduced, dentary with deep horizontal ramus, tall coronoid process and well developed 
angular process, distal humerus articulation less laterally offset, and proximal femur with well-developed trochanters and straight 
shaft. 

Differs additionally from myzopodids in its premaxilla not fused to rostrum and palatine process not well developed, rostrum rela- 
tively wider and with short postorbital process, cranium without basisphenoid pit/s, braincase uninflated, P2 tiny or absent, P3-4 
larger, M1-2 with double ectoflexi, large p3 with two roots (and not p2>p3<p4), p4 with metaconid, m1-3 nyctalodont, talonid rela- 
tively narrower, distal humerus without bilobed medial process (epitrochlea). 

Differs additionally from noctilionoids (noctilionids, mystacinids, mormoopids, phyllostomids, thyropterids, furipterids) in its pre- 
maxilla not fused to rostrum and palatine process not well developed. 

Differs additionally from noctilionids in its hard palate not extending posteriorly beyond level of M3, I2 in toothrow, more than one 
upper premolar, M1-2 without enlarged posteriorly directed heels, mesostyle not protruding buccally, three lower incisors, three 
lower premolars, p4 roots in axis of toothrow, m1-3 nyctalodont, talonid narrower, cristid obliqua does not reach lingual border of 
tooth, coronoid process of dentary taller than condyloid process, ventral margin of horizontal ramus less convex, proximal humerus 
with hemispherical head, distal humerus with less expanded medial process. 
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Differs additionally from mystacinids in its short, wide rostrum, postorbital process, uninflated braincase, two upper (small) inci- 
sors, M1-2 without double ectoflexi, three lower incisors, three lower premolars, m1-3 nyctalodont, dentary with tall ascending 
ramus, distal humerus with spherical rather than narrow central capitulum, broad medial process and short spinous process, prox- 
imal femur with less proximally extended trochanters. 

Differs additionally from mormoopids in its cranium with short postorbital process, uninflated braincase/taller rostrum, M1-2 
without posteriorly enlarged heels, without double ectoflexi, p4 with metaconid, p4>p3>p2, coronoid process of dentary tall, taller 
than condyloid process, m1-3 talonid narrower, proximal femur with well-developed/unreduced trochanters. Differs additionally 
from Pteronotus spp. in m1-3 nyctalodont. 

Differs additionally from phyllostomids in its p4 with tall metaconid. Differs additionally from frugivorous stenodermatines and car- 
olliines, and sanguivorous desmodontines, in its unmodified dilambdodont dentition, and from nectivorous glossophagines, brachy- 
phyllines and phyllonycterines in its short, broad rostrum and unreduced teeth. 

Differs additionally from thyropterids in its cranium with wider rostrum and postorbital process, P2 absent, larger c1, m1-3 talonid 
relatively narrower, and from furipterids in its large C1 and p3 with two roots. Differs from both in its uninflated braincase, dentary with 
tall coronoid process much higher than condyloid process, angular process digit-shaped, and deeper horizontal ramus. 

Differs additionally from extinct speonycterids in its M1-2 not square in occlusal outline, without broad heel, its large p3, m3 with 
hypoconulid, dentary with deep horizontal ramus, and femur with unreduced trochanters. '*° 

Differs additionally from vespertilionoids (vespertilionids, natalids, miniopterids, cistugids, molossids) in its premaxilla not fused to 
rostrum and palatine process not well developed. 

Differs additionally from vespertilionids in p4>p3>p2, proximal humerus with greater tuberosity less proximally extended and less 
well-developed secondary articulation with scapula, distal humerus with spherical capitulum, proximal femur with widely splayed 
trochanters. Differs additionally from the Early Eocene vespertilionids Premonycteris vesper and Sonor handae in its nyctalodont 
m1-3, and additionally from S. handae in its p3 having two roots. 

Differs additionally from natalids in its short, wide rostrum, uninflated braincase, orbit not reduced and infraorbital canal not elon- 
gated, P2 tiny or absent, p4>p3>p2 (not p4=p3=p2), m1-3 talonid narrower, dentary with condyloid process lower than coronoid pro- 
cess and deep horizontal ramus, proximal humerus with less well-developed secondary articulation with scapula, proximal femur not 
laterally bent and trochanters not greatly reduced. 

Differs additionally from miniopterids in its uninflated braincase, deeper, shorter and wider rostrum, 11 and 12 of similar size, C1 
large, M1-2 without posteriorly directed heels, p4 large, p4>p3>p2, m1-3 talonid narrower, deeper horizontal ramus, tall coronoid 
process with condyloid process lower than coronoid process, proximal humerus with less deep fossa for supraglenoid tuberosity 
of the scapula and head hemispherical, distal humerus with spherical capitulum and only short spinous process, proximal femur 
with widely splayed trochanters. 

Differs additionally from cistugids in its large p3 with two roots (not p4>p3<p2) and m1-3 nyctalodont. 

Differs additionally from molossids in its cranium with postorbital process, two (small) upper incisors, three lower premolars with 
roots not oblique in toothrow, tall coronoid process of dentary, only rudimentary second articulation in shoulder joint, distal humerus 
with spherical capitulum and short spinous process, and proximal femur lacking hook-like process on greater trochanter and lesser 
trochanter less proximally extended. 


DESCRIPTION 

Cranium 

General outline and proportions (Figures 1 and S1; Table 1) - The cranium is small (GSL 18.29 mm without premaxilla; Table 1) and 
gracile. It is partly damaged anteriorly (anterior snout broken), ventrally (right pterygoid wing missing) and posteriorly (left petrosal, left 
exoccipital and ventral two thirds of left side of supraoccipital missing). In dorsal view (Figures 1B and S1A), the cranium is slightly 
constricted in the interorbital region, with braincase and rostrum of approximately equal length. The braincase enlarges immediately 
posterior to the orbits for the first half of the braincase, then is rather quadrate in outline. The maximal width is attained in the posterior 
part of the zygomatic arches and, without them, at the posterior extremity of the cranium (across mastoids). Interorbital width is 
approximately one third maximal cranium width (zygomatic), and one half maximal width excluding the zygomatic arches (mastoid). 
In lateral view (Figures 1C and S1B), the cranium is rather flat dorsoventrally; there is no marked dip in height in the interorbital region, 
nor any conspicuous change in basicranial angle, i.e., between palate and basicranium. Except for a slight median bulge, the dorsal 
outline is straight in lateral view and the maximal height of the cranium occurs posteriorly; rostral height is only slightly lower than 
braincase height (Table 1). The occipital surface is flat and forms a marked angle with the dorsal surface. Ventrally (Figures 1A 
and S1C), the petrosal is relatively large (>25% mastoid width) and the promontorium does not well reflect the shape of the cochlear 
canal. 

Rostrum (Figures 1B, 1C, 1E, S1A, S1B, and S1G) - Dorsally, a pentagonal frontal shield is defined by low supraorbital ridges 
ending in postorbital processes and slightly concave frontals with two supraorbital foramina and anteriorly directed grooves that 
issue from these foramina. The point of junction of the ridges is posterior to maximum interorbital constriction (minimum interorbital 
width), the crests formed between them meeting in an angle of ~40 degrees. The anterior frontal fossa at the point of bifurcation of the 
supraorbital ridges appears to be only slightly depressed. Between and anterior to paired supraorbital foramina, is a small single fo- 
ramen. Rostral inflations are present laterally but poorly inflated. The infraorbital canal runs from a large foramen in the anteroventral 
floor of the orbit, exiting anteriorly in an equally large foramen dorsal to P4. At the rostrum’s anterodorsal margin, the nasals do not 
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meet at their most anterior point, giving the anterior margin a scalloped appearance. In anterior view, the nasal aperture is almost as 
tall as wide (Figure S1G). 

Premaxilla (Figures 1F and S1G) - The premaxilla is known from an isolated specimen, UM-VIE-668. The premaxilla body is just 
deeper than the long roots of its two incisors 11 and I2. The nasal process, which appears to have been of constant width throughout 
its depth, is directed posterodorsally but broken caudally with the posterodorsal part missing. Nevertheless, judging from rostrum 
UM-VIE-251, the right and left premaxillae did not meet dorsally, with the nasals dorsally bridging the gap between them. The pre- 
maxilla articulates with the facial portion of the maxilla along a relatively long maxilloincisive margin and, based on UM-VIE-251, prob- 
ably amuch shorter nasoincisive margin. The bone margins of UM-VIE-668 and UM-VIE-251 (and UM-VIE-250) indicate that the max- 
illoincisive articulation was not fused or sutured (no interlocking of the intervening bones '“°), but fibrous/ligamentous. Left and right 
premaxillae appear to have approached or been in contact with each other at the interincisive suture ventrally, and probably loosely 
connected to each other by ligaments. At the point of contact, a short midline flange (which may be broken caudally) is directed to- 
ward the anterior margin of the palate (Figure S1G). This could suggest at least partial ossification of a medial flange of the palatine 
process. However, there is no evidence of a process extending from the anterior margin of the palate and, while some bony division of 
the median incisive fissure might have occurred anteriorly, it evidently did not occur posteriorly. The premaxilla appears to have been 
proclivous, its anterior margin extending anteriorly beyond the level of the upper canine root. 

Palate (Figures 1A and S1C) - The hard palate is relatively short and broad. The tooth rows are roughly parallel anteriorly (C1-P3) 
but from P4 they diverge posteriorly, such that the width of the posterior palate is much greater than the anterior palate. The anterior 
edge of the palate is v-shaped, with the deepest part of the V (its more posterior point) being level with the midpoint of C1 length; there 
is no median spine. The palate’s posterior border has two small lateral indentations that extend anteriorly to the level of the M3 pro- 
tocone. The posterior extension of the midline is level with the posterior face of M3. The maximum length of the palate extends to this 
level. The palate is slightly concave posteriorly (in anteroposterior profile). Palatal foramina occur medial to the anterior face of M1 and 
M2-3. The posterior margin of the palate appears scalloped, with a median projection marking the posterior extent of the vomer. 

Orbit (Figures 1C and S1B) - The orbit is rather large, and bears several openings in its anterior floor. Major ones are identified 
hereafter. There is a concavity in the anteorbital bar, ‘°° of moderate size, here identified as the lacrimal fossa for the lacrimal sac. 
It quickly extends straight through the maxillary bone and opens in the lateral wall of the nasal cavity; this canal is the nasolacrimal 
duct. This configuration is found in other mammals. '°' Ventromedially to the lacrimal fossa is a wide maxillary foramen: it narrows 
anteriorly, then rises a bit dorsally (under the anteorbital bar) and exits anteriorly to the anteorbital bar, through the infraorbital fora- 
men. The latter foramen is therefore at the other side of the lacrimal fossa on the anteorbital bar. On the dorsal surface of this bar are 
also two small pits, which could be the lacrimal foramina. The anteroventral surface of the orbit, posteriorly to the maxillary foramen 
and above the molars, is pierced by several alveolar foramina. There are also three openings medially to the lacrimal fossa and to the 
maxillary foramen. One is large, wider than tall, of irregular outline, and is located just at the transition between the horizontal and 
transverse processes of the palatine (sensu'°’), and just posteriorly to the maxillary. A second opening, small and rounded, is 
retrieved just laterally to the anterior end of the first. Both open posterodorsally. The third, of same shape and size, is posterior to 
the second: it lies just ventrolaterally to the internal choanae and opens posteriorly. The first is the orbitonasal foramen. It is a rather 
uncommon feature in mammals: it is described in insectivorous mammals '°*:'°* where it is a foramen of similar position. According to 
Barone”! and to Giannini et al.,'°* the second opening is the sphenopalatine foramen and the third is the minor palatine foramen 
(posterior palatine foramen in Barghoorn'°°). We thus disagree with some of the foramina Barghoorn (figure 6)'°° recognized: his un- 
named foramen is the sphenopalatine foramen, and his sohenopalatine foramen is rather the lacrimal fossa. The third foramen is the 
posterior opening of a short palatine canal of anteroposterior orientation. This canal opens anteriorly through the major palatine fo- 
ramen between M2 and M3. 

Interorbital and internal nares regions — The narrowest point (Figures 1A, 1B, S1A, and S1C) of the interorbital region is located at 
the level of M3, just anterior to the internal nares. The point where supraorbital crests unite is located posteriorly to this narrowest 
point, at the middle point of the broken zygomatic arch (Figure S1A), and forms the slight bulge of the dorsal outline in lateral 
view (Figure S1B). In lateral view (Figure S1B), the?ethmoidal foramen is marked, just posteriorly and ventrally to the supraorbital 
ridge. The foramen for the anterior opening of the orbitotemporal canal is also present, small, located dorsally to the ethmoidal fo- 
ramen and at the narrowest point of the interorbital region. In ventral view (Figure S1C), the posterior part of the palate is salient 
ventrally and forms a transverse ridge or postpalatine torus between the M3s. Posteriorly, the pterygoids and palatines together 
form raised and straight pterygoid wings, gradually decreasing in height posteriorly, and with a small hamular process. The following 
general observations can be made about the internal nasal skeleton. In the nasal cavity, the ossified nasal septum is well preserved, 
as are the roots of the turbinals on the internal surfaces of the ethmoid, nasal and maxillary bones. Ossified remnants of the fragile 
turbinals themselves and the laminae are also present. Posterodorsally, three to four ethmoturbinals are present in a relatively large 
ethmoturbinal recess. The lamina horizontalis is present and extends laterally, separating the ethmoturbinal recess from the naso- 
pharyngeal passage ventrally. A maxilloturbinal is also apparent, and possibly a frontoturbinal and/or lamina semicircularis. A 
detailed comparative study of Vielasia’s internal nasal skeleton and cribiform plate is being prepared. 

Zygomatic arch (Figures 1A-1C and S1A-S1C) -The entire left zygomatic arch is missing, as is the jugal of the right zygomatic arch. 
In dorsal view, the zygomatic arch diverges posteriorly and seems only curved at its squamosal part, being otherwise rather straight. 
In lateral view, the zygomatic arch originates dorsal to the tooth row and rises posteriorly along a straight line. Its posterior extremity is 
continuous with a ridge running along the squamosal, the mastoid portion of the petrosal and the paroccipital process. 
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Cranial vault (Figures 1B, 1C, S1A, and S1B) - The braincase is rather large but only slightly bulbous laterally. It enlarges quickly 
posteriorly in its anterior half, then only very slightly until its posterior end. There are two large, shallow holes located dorsally, close to 
the sagittal midline, at two thirds of the length. They open in the endocranial cavity on their anterior and lateral sides; these holes may 
be the openings for the transverse sinus (or diploic transverse sinus in Hand'°°). The supraorbital crests are slightly raised, the sagittal 
crest is low to absent anteriorly then taller posteriorly; the lambdoid crest is distinct in its whole length, until reaching the mastoid 
portion of the petrosal. A detailed comparative study of Vielasia’s endocranial cast is in preparation for separate publication. 

Glenoid (Figures 1A and S1C) - The glenoid fossa is flat, and twice as wide as long. A tall postglenoid process forms the posterior 
margin of the articular surface, and it protrudes ventrally, perpendicular to the glenoid fossa. Posterolateral to this process is a large 
postglenoid foramen, the opening of the capsulo-parietal emissary vein. 

Temporal region — Posterior to the glenoid region, the squamosal continues as a curved plate until reaching the mastoid portion of 
the petrosal (Figures S1A and S1B). It shows a ridge, not tall but marked, which is continuous with the lateral edge of the zygomatic 
arch anteriorly and with a shallow ridge on the mastoid portion of the petrosal posteriorly. Its posterior part is the posttympanic pro- 
cess, and it contacts the anterior margin of the paroccipital process of the mastoid portion of the petrosal. Its dorsal suture with the 
parietal is not visible. The parietal contributes to the lambdoid crest but is notched at the lateral extremity of the latter, a part of the 
mastoid portion of the petrosal being therefore visible. The epitympanic recess is present on the medial surface of the posterior 
curved plate of the squamosal. It is located posterodorsally to the glenoid foramen. 

Basicranial foramina, basicranium and occipital region (Figures 1A, 1D, and S1C-S1E) - The anteriormost basicranial foramen is 
the optic foramen. It is large, slightly wider than long and opens in the orbitosphenoid, which constitutes the bony bar separating it 
from the sphenorbital fissure. '°° The latter mostly extends within the alisphenoid. This fissure is ovoid in shape and of constant width. 
The sphenorbital bridge, i.e., the median portion separating the left and right sphenorbital fissures, is very wide and the latter are 
further apart than the two optic foramina. Moreover, the two fissures diverge anteriorly (the sphenorbital bridge widens anteriorly). 
Posteriorly is the oval foramen. It is oval in shape (the major axis of the two foramina converge anteriorly), and of moderate size, 
smaller than the optic foramen. The alisphenoid surrounds it except medially, where a ridge of the presphenoid forms its medial 
wall. This sharp ridge is continuous with the pterygoid wing anteromedially and with a broad and very shallow ridge on the lateral 
aspect of the glenoid fossa posterolaterally. There is no separate round foramen, the latter being most probably confluent with 
the sphenorbital fissure. There are several foramina around the petrosal. One is located anterolaterally to it and is large; it is the pyr- 
iform fenestra. The notch antero-medial to the pyriform fenestra at the anterior end of the transpromontorial sulcus is the carotid fo- 
ramen. The apex of the promontorium contacts the basisphenoid laterally and the basioccipital medially. There is a wide but slim 
aperture between a contact point with the basisphenoid and another with the basioccipital. These apertures are undetermined 
yet. Medially to the petrosal, the long and thin basicochlear fissure separates it from the basioccipital. Finally, there is an almond- 
shaped aperture, located posteriorly to the promontorium (encompassed between the petrosal and the exoccipital); it is the jugular 
foramen. These apertures are not confluent, and clearly separated by connections between the petrosal and other basicranial bones. 
The hypoglossal foramen opens in the exoccipital bone, anteroventrally, and is located just posteromedially to the jugular foramen. Of 
the foramen magnum, only remains its ventral and right margins. It seems ovoid in shape, with a flat ventral base, and is wider than 
high (Figures S1C and S1E). The exoccipital participates to the paracondylar process together with the posterior part of mastoid of 
the petrosal. The occipital condyle is narrow, weakly projected ventrally with a small dorsal expansion. Its external margin is thick and 
shows two concave surfaces, one posteriorly to the hypoglossal foramen, and another postero-medially to the paracondylar process. 
There is a groove between the posterior aspect of the petrosal and the ventral margin of the foramen magnum, throughout the ba- 
sioccipital and the exoccipital. This groove is deeper between the occipital condyle and the paroccipital process, at the condyloid 
fossa. Dorsally, the supraoccipital forms marked and developed lambdoid crests when contacting the parietals. These crests, in oc- 
cipital view (Figure S1E), are straight and diagonal in orientation; the lateral bulge of the cranial cavity is visible in this view. In lateral 
view, the supraoccipital is only slightly concave in its dorsal half, giving a general straight and vertical shape to the posterior face of 
the cranium. There is a large aperture in its most lateral part, where lies the mastoid portion of petrosal. This aperture follows the 
anterior and medial edges of the mastoid but do not contact it, there is a wide gap in the location of the mastoid foramen, this struc- 
ture could be refered to as the “mastoid aperture.” 

Petrosal (Figures 1D and S1B-S1E; Data S1F) - The auditory bulla, stylohyal, and auditory ossicles are missing. Only the right 
petrosal is preserved in situ (Figure S1C). Its pars cochlearis contacts the basicranium through four bone contacts (see above). 
On the ventrolateral surface, the fenestra cochleae and fenestra vestibuli (ventromedially and dorsolaterally respectively) are large 
openings (Figure S1D). The fenestra vestibuli is oval, longer than wide. The fenestra cochleae is bean-shaped, being much wider 
than long. Both are separated by a large crista interfenestralis. Just anteriorly to the fenestra vestibuli is an aperture. In Chiroptera, 
except pteropodid bats, ‘°° the facial nerve does not run in a canal then a sulcus as in some other mammals (e.g., marsupials, '°’ 
artiodactylans, ‘°° notoungulates '°°). This aperture therefore corresponds to both the hiatus Fallopii (the exit foramen for the greater 
petrosal nerve branch of facial nerve VII) and the secondary facial foramen (where, typically, the rest of the facial nerve VII exit and 
continues along in the facial sulcus). This foramen is hereafter named a facial foramen. The facial sulcus then widens distally: it is 
firstly a short and narrow groove then a wider area. The posterior end of this area, just dorsolaterally to the base of the fenestra ves- 
tibuli, is a short and thick rim. This is a vertical thickening of the facial sulcus floor, and it is pierced by asmall foramen which opens ina 
sinus of the squamosal. A large concave surface follows this rim posteriorly and ventromedially, posterior to the windows; it is very 
wide, with a bean shape. In this surface are a large and excavated stylomastoid notch dorsolaterally and a shallow stapedial muscle 
fossa ventromedially. The previous short and thick rim is continuous anteriorly with the floor of the fenestra vestibuli and posteriorly 
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with a bony expansion. The latter corresponds to the tympanohyal attachment site. It is directed antero-ventrally and its apex is large 
and concave. A very subtle transpromontorial sulcus runs roughly parallel to the facial sulcus on the whole length of the promonto- 
rium, at the antero-posterior level of the fenestra vestibuli. There also is a very shallow sulcus for the stapedial artery, just medially to 
the fenestra vestibuli, joining the transpromontorial sulcus at a right angle. Posteriorly, the transpromontorial sulcus passes just pos- 
teriorly to a tall bony spike. Anteriorly to this spike is a small crest that runs along the ventromedial edge of the petrosal and joins 
the most medial bony contact with the basioccipital. A similar spike or tubercle occurs in natalids (see figure 5 of Morgan and Cza- 
plewski 18°), emballonurids (JM pers. obs.) and molossids (see figure 12 of Henson et al.*°), and also Onychonycteris finneyi (see sup- 
plemental figure 2 of Simmons et al.°) In natalids, it bears a small hook-like process that in life partly encircles the internal carotid 
artery and contacts part of the posterior end of the stylohyal.'©° In emballonurids and molossids, it is one of the attachment sites 
of the auditory bulla and it greatly resembles that described here; we propose a similar function. The shape, location and putative 
function of this spike are very similar to those of the metatherian rostral tympanic process (as illustrated by Sanchez-Villagra and 
Wible'°’). Being regarded as a typical marsupial feature, we do not name it in the same way, but this should be emphasized. 

There is no epitympanic wing. The antero-lateral surface of the promontorium shows a slight depression corresponding to the area 
(a “fossa” when more marked) for the tensor tympani muscle insertion. Ventromedially to its posterior ridge is the anterior end of the 
transpromontorial sulcus. The bulge of the cochlear canal is decipherable on the promontorium surface, but is not very well outlined; 
there are at least two turns in the cochlear canal, with two and one-quarter turns according to Davies et al.'° or two and a third turns 
according to Ekdale.®° There is a small caudal tympanic process, tubercle-shaped, immediately posteriorly to the fenestra cochleae. 
Posteriorly to the fenestra cochleae and the caudal tympanic process is a concavity, which corresponds to the petrosal contribution 
to the anterolateral margin of the jugular foramen. There is a great bony expansion dorsolaterally to the floor of the facial sulcus, which 
runs anteriorly with a spike directed anteroventrally; it is the tegmen tympani. It is pierced by a foramen of comparable size to that of 
the fenestrae and which opens posteroventrally; it is probably the foramen for the superior ramus of stapedial artery. '°' The posterior 
base of the tegmen tympani contributes to the tympanohyal attachment site. The epitympanic recess is not located on the petrosal 
but on the squamosal (see above), just dorsolaterally to the tegmen tympani. It is quite deep. There is a small but definite fossa in- 
cudis, at the posterior extremity of the dorsolateral edge of the tegmen tympani. This fossa is oriented toward the epitympanic recess 
of the squamosal, ventrolaterally. The mastoid portion of the petrosal bone is massive, of comparable size to that of the pars coch- 
learis. It is widely exposed ventrally, occipitally, and laterally. The lateral edge of the mastoid forms a large and blunt ridge, the pa- 
roccipital process, continuous with the ridge of the squamosal anteriorly (the posttympanic process) and the paracondylar process 
posteriorly (Figures S1B, S1C, and S1E). The dorsal part is widely visible occipitally and less so laterally, within the “mastoid aper- 
ture” (Figure S1E). 

Bony labyrinth (Figures 3A and S1F; Data S1F) - The cochlea has a volume of 3.05 mm? and comprises ~74% of total inner ear 
volume. The cochlear canal makes 2.25 turns according to the criterion of Davies et al.,'° or 2.3 turns following Ekdale.®° It has a 
maximum width of 2.50 mm (first turn diameter). The aspect ratio of the cochlear spiral in profile®® is 0.62; the fenestra vestibuli 
(oval window) is elliptical (stapedial ratio 1.75). A secondary bony lamina is present along the radial wall of the cochlea canal and 
appears to extend well beyond the basal turn, almost reaching the apex. The length of the cochlea canal/basilar membrane is 
~12.38 mm (method of Davies et al.'°). The plane of the basal turn of the cochlea forms an angle with the plane of the lateral semi- 
circular canal of ~25 degrees. The canaliculus cochleae for the cochlear aqueduct is small and slightly curved. The dorsal edge of the 
cochlear fossula bears a horseshoe shaped keel extending anteriorly on the dorsal surface of the cochlear canal and connecting 
medially to the root of the cochlear aqueduct. The spherical recess is not distinguishable from the elliptical recess in the vestibule. 
Vielasia possesses an anterior excavation for the anterior and lateral ampullae, expressed as a pedestal for the anterior and lateral 
ampullae on the endocast, as well as a posterior excavation for the posterior ampulla and common crus. The common crus is tall and 
slender but thicker than the anterior, posterior and lateral canals. The bony channel for the vestibular aqueduct leaves the inner ear 
directly medial to the vestibular aperture of the common crus. The posterior limb of the lateral semicircular canal opens into the pos- 
terior ampulla. The angle between the planes of the anterior and posterior semicircular canals (101 degrees) is greater than between 
the lateral and posterior canals (89 degrees) which is similar to that between the anterior and lateral canals (90 degrees). The radius of 
the arc (R'2“) of the anterior semicircular canal (0.92) is greater than the posterior canal (0.77) which is greater than the lateral semi- 
circular canal arc (0.74). 

Hearing frequency range reconstruction (Figure S4; Data S1K and S1L) - Davies et al. “` performed several multiple regression an- 
alyses to assess the relationship between echolocation call parameters (minimum, peak, and maximum frequencies) and putative 
hearing-related measurements (basilar membrane length, number of cochlear turns, body mass). Their best-fitting model involved 
only cochlear measurements (basilar membrane length, number of cochlear turns), with empirical formulae provided”? (see also addi- 
tional file 1, table s6 of Davies et al.'°). Using these formulae, we reconstructed echolocation call parameters for Vielasia sigei, given 
its basilar membrane length and number of cochlear turns (Data S1F and S1G), as minimum echolocation frequency of 30.19 kHz, 
peak echolocation frequency of 41.95 kHz, and maximum echolocation frequency of 55.57 kHz. Davies et al.'° also inferred low and 
high limits of hearing frequencies at 30 dB and at 60 dB, using the methods of West'°* and Manoussaki et al.'°° They performed OLS 
regressions of (1) logi9 low limit of hearing frequencies against log;9 of the product of basilar membrane length and number of 
cochlear turns, and of (2) log; high limit of hearing frequencies against log;9 of relative basilar membrane length (i.e., ratio between 
basilar membrane length and cube root of body mass). Davies et al.'° used a dataset containing 38 mammal species but discarded 
two subterranean rodents and two marsupials for performing OLS regressions. The latter dataset contains 14 chiropteran species 
and 20 non-bat placental mammal species. Davies et al.'° performed OLS regressions both with and without the 14 bat species, 
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and therefore obtained two OLS regression lines for each of the four cases (low/high limit, 30/60 dB). We re-analysed the data (Data 
S1K; Figure S4) and estimated low and high limits of hearing frequencies at 30 dB and 60 dB for Vielasia sigei, comparing these also to 
values calculated by Davies et al.'° for the theoretical “chiropteran ancestor” (Data S1L). 

Eye size — Eyeball size in Vielasia sigei was estimated as ~2.3 mm diameter or less, using the method of Brooke et al.°° (and see 
Thiagavel et al.?? for validation in bats). Using this method, a virtual spherical eye ball was fitted into the socket in a model of the ho- 
lotype UM-VIE-250 without touching any part of the skull while the gap between bone and eyeball on the median side was maintained 
at ~1 mm. 

Upper dentition 
Dental terminology follows Legendre (see figure 1),'°* Hand, '°° and see Figures S1H and S11. 

Upper incisors - 11 and I2 are known from isolated teeth and are present in premaxilla UM-VIE-668 (Figures 1F and S1G). Both I1 
and I2 are small, of similar size, hypsodont, bilobed with the anterior lobe longer and taller than the posterior lobe, and lacking cingula. 
The incisors in UM-VIE-668 are worn, but 11 has lobes of more equal size (both length and height) than 12 in which the posterior lobe is 
significantly smaller such that the tooth is nearly caniniform. The root in both 11 and 12 is long, at least 2.5 times crown height. Upper 
incisors are not particularly reduced, with 11 length in toothrow approximately half M3 length. 

Upper canine - In cross-section C1 is anteroposteriorly elongate (Figures 1A, S1C, and S1H). It has a tall, anteriorly placed single 
cusp (paracone) with two main crests separating its two faces; the posterior crest is sharper that the anterior crest. The lingual face is 
flat and slightly depressed; the buccal face is convex in section but its distal surface is somewhat flatter. A discontinuous cingulum is 
well developed around most of the crown but in some specimens attenuates anterobuccally and is absent anteriorly. A small but 
distinct cingular cuspule is formed where the lingual and buccal cingula meet the cusp’s posterior crest. 

Upper premolars - UM-VIE-250 and UM-VIE-251 preserve alveoli for two premolars (Figures 1A, S1C, and S1H), interpreted here 
as P3 and P4. The first upper premolar (P2) is absent in UM-VIE-250 and UM-VIE-251 but an alveolus for a very small, single-rooted 
tooth in this position is present in UM-VIE-252. P3 is known from several specimens (Table 1). It is a relatively small, three-rooted 
tooth, round to triangular in occlusal shape. It has a single central cusp that is rounded in section anteriorly and has a sharp posterior 
crest. A cingulum is well developed around the crown but dissipates buccally where the central cusp is most convex/at the cusp’s 
greatest convexity. P4 is a three-rooted tooth that is conspicuously larger than P3. It is slightly wider than long and T-shaped in 
occlusal view. It has a tall and sharp buccal cusp (paracone), a low and weakly cuspate parastyle that projects slightly anteriorly 
and alow rounded and expanded lingual lobe or heel. It has a distinct cingulum along the buccal, anterior, lingual and posterior sides; 
the anterior cingulum dissipates/is absent in its medial part. The lingual cingulum increases in height anteriorly where it bears a 
distinct lingual cuspule (protocone). A curved posterobuccal crest (postparacrista) runs from the paracone to the posterobuccal 
corner of the crown. 

Upper molars - There are three upper molars (Figures 1A, S1C, and S1H). M1 is a three-rooted tooth that is much wider than long 
(mean ratio 1.4; Table 1). The metacone is taller than the protocone which is taller than the paracone. The buccal border has an overall 
oblique direction with respect to the paracone, and has a single, long, shallow pre-ectoflexus anterior to the mesostyle. The well- 
developed, crested parastyle projects anterolingually. The postpara- and premetacrista meet close to the buccal cingulum, with 
the mesostyle no more than a slight elevation of the buccal border crest. The postparacrista and premetacrista meet at an angle 
of ~60°. The postmetacrista is longer than the premetacrista which is longer than the subequal pre- and postparacrista. The proto- 
fossa is closed posteriorly by the postprotocrista which meets the base of the metacone, or in some specimens also extends beyond 
the metacone and connects with the posterior cingulum. The preprotocrista extends along the anterior flank of the paracone to meet 
the base of the parastyle in an acute angle. A paraloph, extending from the base of the paracone to the preprotocrista, is variably 
present. There is no metaloph. The protocone is relatively tall, narrow (anteroposteriorly), and lingually extended. The basal posterior 
face of the protocone extends lingually and slightly posteriorly, developing as a rounded heel which has a distinct but weak cingulum, 
best developed anteriorly and posteriorly. There is no hypocone, but a swelling in the posterolingual cingulum is sometimes present 
(e.g., UM-VIE-250, Figures 1A and S1C). In UM-VIE-250 a depression in the heel’s posterior cingulum appears to represent a wear 
facet for the protoconid of m2. M2 has a similar structure to M1 but differs as follows. M2 is wider than M1 (and is proportionately 
wider; mean width to length ratio 1.61) and its protocone heel shorter anteroposteriorly. The buccal border, from anterior to posterior, 
is directed anteroposteriorly rather than obliquely, and this border also has a deeper and longer ectoflexus, more medial in position 
than in M1. The paracone and metacone are more equal in size, their buccal valleys more similar in area, and the preparacrista, post- 
paracrista and premetacrista more equal in length although still shorter than the postmetacrista. The protocone heel is less poste- 
riorly and more lingually extended, and the basal cingulum more or less continuous around the protocone. M3 is conspicuously wider 
than long, but narrower than M2 and typically subequal to M1. It retains a metacone, which is lower in height than the paracone and 
protocone. The preparacrista is longer than the postparacrista which is just longer than the premetacrista; there is no postmetacrista. 
The tooth has a well-developed, cuspidate and anteriorly directed parastyle. The paracingulum is narrow, as in M1-2, and the pro- 
tofossa deep. The posterobuccally directed postprotocrista meets the base of the metacone, closing the protofossa. There is no 
cingulum around the base of the protocone, and no buccal, anterolingual or posterior cingulae. 

Dentary and lower dentition 

Dentary - The horizontal (mandibular) ramus (UM-VIE-254 and UM-VIE-255) has a slightly convex ventral margin, best seen in UM- 
VIE-254 (Figures 1G-1H and S11), but is of uniform medium depth (2.1 mm), approximately equivalent to 1.5 times m1 height. Its sym- 
physis extends posteriorly to mid-p3 (Figures 11 and S11). The anterior mental foramen occurs below the incisors; the posterior mental 
foramen occurs below c1-p2. The dentary has a tall ascending ramus, with the condyloid process occurring at a level well dorsal to 
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the tooth row but conspicuously lower than the tall coronoid process. The anterior edge of the coronoid rises nearly vertically, its tip is 
missing but is at least twice tooth crown height. The well-defined angular process is digit-shaped, more or less in line with the hor- 
izontal ramus and occurs at the level of the tooth alveoli. An elongate mandibular foramen opens well ventral to the tooth row. The 
dental formula is i3, c1, p2, p3, p4, m1-3, where p2 is single rooted and p3-m3 have two roots each. 

Lower incisors — No lower incisors are preserved in situ but alveoli in UM-VIE-254 and CT data indicate that three incisors were 
present (Figures 1G-11 and S11). Isolated lower incisors (e.g., Figure S11) are small tri-lobed teeth, longer (in toothrow) than wide, 
with wide, well-developed basin, brachyodont, without cingulid. 

Lower canine — The c1 is anteroposteriorly elongate, and there is a single, tall primary cusp (protoconid) and posteriorly-extended 
heel (talonid) with a small transverse basin (Figures 11 and S11). The protoconid’s buccal face is convex, while the posterior face is flat 
or slightly concave. A cingulid is very well developed lingually but is indistinct or absent buccally. Anterolingually, the cingulid joins the 
vertical anterior crest of the principal cusp in a sweeping curve that extends up the anterolingual face of the canine attenuating at 
about half cusp height. Posterolingually, the cingulid is lower but still pronounced. A small cingular cuspid is present posteriorly. 
The root is of similar proportions to the crown. 

Lower premolars — The first lower premolar p2 is a single rooted tooth (Figures 11 and S11). It has a low main cuspid, is longer than 
wide, elliptical in occlusal view and triangular in lateral view. Cristids run anteriorly and posteriorly from the tip of the main cuspid 
(protoconid) to small cingular cuspids at the crown’s most anterior and posterior points, respectively. The lingual cingulid is stronger 
than the buccal cingulid but both are distinct. The middle lower premolar p3 is a larger two-rooted tooth (Figures 1G-11 and S11). In 
occlusal shape and general proportions, the crown is rectangular, and longer than wide. It consists of a tall main cuspid (protoconid). 
This cuspid has a cristid running from its tip anterobuccally to a well differentiated paraconid, and a less distinct cristid running pos- 
teriorly to a cingular basin. No metaconid is present. There is a discontinuous basal cingulid around the crown which is better devel- 
oped lingually than buccally where it attenuates approximately midway along the length of the crown, at the greatest convexity of the 
protoconid. A small but distinct cingular cuspid occurs at the most posterior point of the crown. The posterior lower premolar p4 is a 
two-rooted tooth similar in length to p3 but conspicuously wider. It consists of a tall main cuspid (protoconid) with a cristid running 
from its tip anterobuccally to a well differentiated paraconid and a posterior cristid to a less distinct metaconid, the former being 
smaller and closer than the latter to the main cuspid (protoconid). Variably, a posterior (entocristid-like) cristid continues from the 
metaconid to the posterolingual cingulid. A basal cingulid is more or less continuous around the crown. A small cingular cuspid pro- 
jects posteriorly. 

Lower molars — In m1, the trigonid is shorter but slightly wider than the talonid (Table 1). The trigonid basin is significantly higher 
than the talonid basin (Figures 1G, 11, and S11). The protoconid is the tallest cuspid followed closely by the metaconid which is taller 
than the subequal paraconid, hypoconid and entoconid. The hypoconulid is distinct but low; it is lingually situated and the postcristid 
runs from the hypoconid directly to the hypoconulid, in the nyctalodont pattern.'©° In occlusal view (Figure 1H), the metacristid is 
straight and almost at right angles to the long axis of the tooth, whereas the paracristid is longer and slightly more curved, reflecting 
a slightly anteriorly directed paraconid. In posterior view, the angle formed by the metacristid components is less than 90°, whereas 
that formed along the paracristid is conspicuously more than 90°. The angle formed between the paracristid and metacristid at the 
protoconid is approximately 50°. The cristid obliqua joins the posterior wall of the trigonid (postvallid) at the centre of the crown, 
below the junction of the two components of the metacristid, such that the hypoflexid is relatively deep. The anterior, buccal, and 
posterior cingulid is continuous, but is thinnest around the hypoconid. There is no lingual cingulid. The m2 is described in so far 
as it differs from m1. The trigonid is more anteriorly compressed such that the angle formed between the paracristid and metacristid 
at the protoconid is approximately 40°. The lengths of the paracristid and metacristid are more equal than in m1. The talonid cuspids 
are relatively shorter, such that there is a more conspicuous difference in height between trigonid and talonid cuspids in m2. 
Compared to m1-2, m3 is a narrower tooth (Table 1), with talonid width conspicuously less than trigonid width. The talonid cuspids 
are shorter, such that there is more difference in height between trigonid and talonid cuspids than in m1-2. The angle formed between 
the paracristid and metacristid at the protoconid is approximately 35°. The basal cingulid is briefly absent around the greatest con- 
vexity of the hypoconid. 

Postcranium 
Terminology and orientation for skeletal element morphology follows Vaughan?” and Smith.°° 

Scapula - UM-VIE-671 represents a fragment of a left scapula (Figure 1K). Few structures are preserved but those include the 
glenoid fossa, the supraglenoid tuberosity and the bases of the coronoid process (ventral surface) and acromion process (dorsal sur- 
face). The glenoid fossa of the scapula is pear-shaped, being narrower cranially than caudally. It is defined by a low lip except dor- 
somedially where a smooth, poorly developed facet (dorsal articular facet) marks an articulation with the greater tuberosity of the 
humerus. 

Humerus — UM-VIE-672 is the proximal end of a right humerus (Figures 1L and 1M). It is 4.6 mm wide dorsoventrally and 3.5 mm 
deep mediolaterally without pectoral ridge and ~5 mm deep mediolaterally with pectoral ridge; this fragment is 10.7 mm long. The 
greater tuberosity is knob-like and extends proximally just beyond the head. The lesser tuberosity is more wing-like, being markedly 
ventrally extended, and does not extend proximally beyond the head. The head is hemispherical and centred between the tuberos- 
ities; it extends caudad only just beyond the shaft. The head is separated from both tuberosities by equally deep grooves, the groove 
between the head and greater tuberosity being narrower. A shallow fossa is developed anterior to the head and between the two 
tuberosities but closer to the greater tuberosity. The pectoral ridge is long and tall, its edge turned so that its posterior surface is 
concave and its medial surface forms the lateral part of the broad, deep bicipital groove. The pectoral ridge meets the shaft distally 
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~8 mm from the most proximal point of the greater tuberosity. A ventral crest (medial ridge) that runs distally from the base of the 
lesser tuberosity for ~3 mm is anteriorly concave. A less pronounced crest extends from the distal part of the greater tuberosity 
before dissipating. 

The distal end of the shaft is oval in cross-section (~1.8 x 1.3 mm). UM-VIE-673 represents the distal portion of a left humerus (Fig- 
ure 1N). It is 4.4 mm wide mediolaterally and 2.4 mm deep anteroposteriorly; this fragment is 6.7 mm long. The articular surface is 
slightly offset with respect to the shaft so that, in anterior view, both the trochlea and lateral epicondyle rim are just lateral to its medial 
and lateral edges respectively. The central capitulum is spherical and occupies ~40% the articular surface width. The medial process 
(epitrochlea) is relatively broad, representing ~30% the width of the articular facets (including distal spinous process representing 
one-half that width ~17%). The distal spinous process is well separated from the trochlea and short (does not reach the distal artic- 
ular surface). The trochlea and central capitulum are more prominent than the lateral epicondyle. The trochlea extends slightly further 
distally than the capitulum; they are separated by a trochlear groove. A ridge extends from this margin proximally along the shaft to 
enclose, medially, a deep, broad, central radial fossa. 

Radius - UM-VIE-674 represents the proximal portion of a left radius (Figure 10). It is ~2.8 mm wide dorsoventrally and ~2.4 mm 
deep anteroposteriorly; this fragment is 8.2 mm long. The epiphysis is ~1.7 times as thick as the shaft, and inclined craniad (tilting 
from the lateral to medial sides) with respect to the long axis of the shaft. In anterior view, the epiphysis is a rounded triangular shape, 
with the most acute angle being blunt and proximally directed. The surface is concave with the articular facet for the capitulum being 
large and round, occupying ~60% the width of the articular surface. There is no conspicuous ridge separating articular facets for the 
capitulum and trochlea. Medial and lateral winged articular processes are similar in size and distally tilted. The anterior flange of the 
radius articulates with the anterior face of the trochlea, and the posterior flange articulates with the epicondyle anterior face. Poster- 
omedially, there is no clear facet indicating the contact area for the posterior articular facet of the cubitus of the ulna. Immediately 
distal to the anteromedial part of the distal rim of the articular surface is a deep slit (flexor fossa) into which the tendons for the 
Mm. triceps and brachii and brachialis would insert. The distal end of the shaft is circular in cross-section (diameter ~1.65 mm). 

Femur — UM-VIE-675 represents the proximal portion of a left femur (Figure 1P). It is ~3.2 mm wide mediolaterally, ~1.5 mm deep 
anteroposteriorly; this fragment is ~4.8 mm long. It is relatively large, its width across the trochanters being 70% of that across the 
proximal humerus tuberosities. The head of the femur, clearly distinct from the rest of the epiphysis, is spherical, slightly anterome- 
dially directed and has a distinct but shallow fovea capitis. The head is relatively large (diameter ~1.4 mm), slightly wider than the 
distal shaft, and with its medial edge aligned with the shaft’s medial edge. Two well-developed trochanters are located opposite 
each other but both on the posterior surface of femur. The head lies between the trochanters, along the longitudinal axis of the shaft 
and extends anteromedially. The trochanters are similar in size and widely spread. They are connected to one another by a weak 
horizontal ridge across the posterior surface of the femur just distal to the head. The lesser trochanter does not extend proximally 
beyond the base of the head; the greater trochanter extends further proximally (half head height). Grooves separate the head 
from the shaft and the greater trochanter. The shaft appears to have been straight, as least proximally. It is slightly flattened such 
that low crests extend from greater and lesser trochanters. The distal end of the shaft is oval in cross-section (~1.10 x 0.89 mm). 

Comparisons of postcranial elements with other stem bats — Vielasia’s skeletal elements exhibit specializations for flight found in all 
bats, extant and extinct (e.g., reorientation of scapular spine, absence of olecranon fossa, supinator ridge and entepicondylar fora- 
men in humerus, elongated radius with respect to humerus, reorientated femur with poorly-developed neck, prominent and proximal 
head, greater and lesser trochanters equally developed, lack of digital fossa).°4°"*° 

Vielasia’s shoulder joint (with hemispherical head on proximal humerus and scapula with correspondingly simple, pear-shaped gle- 
noid cavity) is of a form regarded as the least specialized type among extant bats, as found in (e.g.) pteropodids and rhinopomatids.°” 
In Vielasia, however, the major tuberosity of the proximal humerus extends proximally beyond the head and a shallow pit receives the 
greater tuberosity of the scapula’s glenoid fossa at stages in wingbeat cycle. These represent specializations that reduce the pro- 
natory movements of the abducted forearm during the downstroke of the wing-beat cycle,°’ possibly facilitating longer time on 
the wing. Among stem bats a second articulation in the shoulder is better developed in Hassianycteris and Palaeochiropteryx spp. 
but not in Onychonycteris and Icaronycteris°*° and possibly not in Archaeonycteris (see figure 23 of Russell and Sigé”°). 

Similarly, the elbow joint morphology in Vielasia is similar to that found in stem bats. Vielasia retains a spherical central capitulum 
and broad medial epicondyle with low spinous process (indicative of a large medial muscle mass), a morphology also seen in Icar- 
onycteris,“° Onychonycteris and Archaeonycteris spp. (see figure 24 of Russell and Sigé’°) and to a less degree in Palaeochiropteryx 
and Hassianycteris spp. (in which the spinous process extends to the level of the distal trochlea; see figure 11 of Russell and Sigé”°). 
This arrangement enables considerable rotational (anteroposterior and lateral) movement in the joint, facilitating the maneuvrable 
flight used today by (e.g.) phyllostomids, hipposiderids and pteropodids. In more archaic bats, the articular surface is laterally posi- 
tioned, with the capitulum either completely or mostly lateral to the lateral edge of the humeral shaft, and the medial edge of the 
trochlea well lateral to the medial edge of the humeral shaft*®; this condition is also seen among extant (e.g.) pteropodids, rhinolo- 
phoids, phyllostomids, some mormoopids, and in noctilionids. In Vielasia, it is a little more specialized: the capitulum lies mostly 
medial to the lateral edge of the humeral shaft, such that the articular surface is displaced only slightly laterally with respect to the 
midline of the shaft. 

With well developed, broad trochanters, low lesser trochanter and slender, straight shaft, Vielasia’s proximal femur is also similar to 
that of other stem bats including /caronycteris and Archaeonycteris spp. (SJH pers. obs.) but differs from the condition in Palaeochir- 
opteryx (see figure 15 of Russell and Sigé”°) in which the head is more proximally directed and the trochanters less splayed. It is 
similar to that of most extant bat species with generalized roosting habits (i.e., hang from branches or ceilings, cling to vertical 
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surfaces and/or rest in crevices) such as pteropodids, vespertilionids and emballonurids.°°:'°? It lacks specializations seen in bats 
that hang pendant-like, such as mormoopids and rhinolophids (trochanters greatly reduced, head directed proximally), or bats 
that walk, such as vampires, mystacinids and some molossids (robust shaft, well developed lateral and medial ridges, proximally 
extended lesser trochanter '°’). 

Other specimens referred to Vielasia sigei - More than 400 craniodental and postcranial fossil specimens recovered from Vielase 
locality in the Quercy Phosphorites, southwestern France°”’°° °° appear to be referable to Vielasia sigei. Qualitative examination, uni- 
variate statistics, mixture analysis and principal components analysis of craniodental specimens in the fossil sample are consistent 
with the material representing a single species - the common Vielase bat species noted by Legendre et al.°° and described here as 
Vielasia sigei n. gen. et sp. A second, much smaller bat species is represented by rare teeth, also noted by Legendre et al.°°; these 
teeth will be described separately and are not included in the sample analysed here. Univariate and multivariate analyses of 131 com- 
plete/measurable craniodental specimens were performed in the statistical program PAlaeontological STatistics (PAST '''°°), with 
univariate summary statistics given in Data S1A. Coefficient of variation (CV) values lie within the range expected from a single pop- 
ulation, i.e., CV values between 4 to 10,'°° and PCA analysis (using iterative inputation for missing data) also found no evidence of 
clustering that would suggest the presence of more than one population/sexual dimorphism among isolated teeth referred to this 
species (Data S1A). For most variables, mixture analysis '’°'’' showed no evidence of clustering or dimorphism in the sample of iso- 
lated teeth. However, for lower first molars (N=41), mixture analysis did detect two groupings in trigonid widths (m1WT) with means at 
1.30 mm and 1.47 mm (prob. 0.536 and 0.464, respectively), and two groupings in talonid widths (m1 Wt) with means at 1.20 mm and 
1.33 mm (prob. 0.541 and 0.459, respectively) (Data S1A). The well-preserved Vielase fossil tooth sample provides an opportunity to 
further examine dental variation, including sexual dimorphism, in an early Eocene bat and is the subject of a separate analysis by 
additional authors. 

Unfused premaxilla and familial distinction 

In Vielasia, the premaxilla contains two relatively large incisors (Figures 1F and S1G), the maximum number found in any bat and the 
number found in all stem bats. However, unlike all other early Eocene bat taxa in which the morphology is thought to be clear, 
Vielasia’s premaxilla was only loosely attached to the rostrum (see Description above), as it is in many extant bats, including in 
both yangochiropterans and yinpterochiropterans. '*° In Vielasia, the maxilloincisive articulation in the adult bat appears to have 
been ligamentous, rather than fused or sutured; left and right premaxilla bodies approach one another anteroventrally, or meet in 
a loose midline interincisive contact; a short medial flange extends caudally from the premaxilla body but a palatine process from 
the anterior margin of palate is absent (we interpret/score this as palatine process of the premaxilla being reduced or absent); nasal 
processes of right and left premaxillae do not meet dorsally (nasals bridging the gap); and the anterior edge of the premaxilla extends 
beyond the level of the upper canine root (proclivous condition’). In bats, the variable form of the premaxilla is taxonomically diag- 
nostic among extant and extinct families, '“° and V. sigei exhibits a combination of premaxillary features that suggests familial distinc- 
tion for Vielasia may be warranted, but a re-examination of this complex in other stem bats is needed before this can be confirmed as 
unique. The adaptive or functional significance of a reduced, relatively mobile premaxilla in bats is currently not understood, but pos- 
sibilities include increased gape, inflight prey capture and manipulation, modulation of echolocation call, noseleaf function, reduced 
interference with returning echolocation call, and/or olfactory function (as noted by Orr et al.'’° and references therein). 

Scoring of premolar positions 

Extant and fossil bats have a maximum of three premolars, which have been traditionally identified as P2-4/p2-4 in descriptions of bat 
comparative anatomy (e.g., Smith et al., Hand et al.,°' Czaplewski et al.,°° Simmons and Geisler,*° Sigé, '°' Russell et al., and 
Hand et al.'““). To facilitate comparisons between these studies and the description of Vielasia sigei presented here, we also refer to 
bat premolars as P2-4/p2-4 in the descriptive and comparative aspects of our study. However, our craniodental character matrix is 
based on that of O’Leary et al.,°' who identified the premolars of bats as comprising P1/p1 P4/p4 and P5/p5 of the ancestral eutherian 
complement of five premolars; to ensure that we have correctly compared equivalent tooth positions between bats and the non-bat 
eutherians in our matrix, we have followed O’Leary et al.’s°' scheme for bat premolars when scoring our dental characters (Morpho- 
bank Project 4407, http://morphobank.org/permalink/?P4407). 

All of the five known maxillary specimens of V. sigei contain two upper premolars (P4 and P5 of O’Leary et al.°' In four of these five 
specimens (including the holotype UM-VIE-250; Figures 1A and S1C), there is no sign of, nor space for, P2 (P1 of O’Leary et al.).°' A 
fifth maxilla preserves a possible alveolus buccal to the centre of P3 (P4 of O’Leary et al.°') which could represent either a tiny, 
buccally extruded, variably present premolar or a dental atavism (Figure S1H). However, in our phylogenetic analysis, we scored 
the anteriormost premolar as being absent in Vielasia sigei, as this is clearly the modal condition. 


PHYLOGENETIC ANALYSES 

In order to test the evolutionary relationship of Vielasia sigei to other chiropterans (extant and fossil), we carried out a series of phylo- 
genetic analyses using both morphological and total evidence (= combined morphological and molecular) matrices, as follows. 
Morphological matrix 

Our morphological character list comprised 2665 craniodental characters, which is a reduced version of the 4541 phenomic char- 
acter matrix of O'Leary et al.°' (with postcranial, soft tissue, behavioural and developmental characters deleted). We used this matrix 
as a starting point because it is the largest (in terms of characters) morphological matrix currently available for resolving higher-level 
relationships within mammals, and because it has already been scored for the early Eocene bats Onychonycteris finneyi and Icaro- 
nycteris index, plus six extant bats that collectively represent all major crown chiropteran lineages, namely: the yinpterochiropteran 
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Pteropodidae (Pteropus giganteus) and Rhinolophoidea (the rhinopomatid Rhinopoma hardwickii), and the yangochiropteran Embal- 
lonuroidea (the emballonurid Saccopteryx bilineata and the nycterid Nycteris thebaica), Vespertilionoidea (the vespertilionid Myotis 
lucifugus), and Noctilionoidea (the mormoopid Pternotus parnelli) O’Leary et al.°' As outgroups, we retained from O'Leary et al.’s°' 
original matrix the stem eutherians Maelestes gobiensis, Ukhaatherium nessovi, and Zalambdalestes lechei, plus five plesiomorphic 
members of the same placental superorder to which Chiroptera belongs, namely Laurasiatheria: the early Paleocene Protungulatum 
donnae (a probable stem-euungulate), the early Eocene carnivoramorphan Vulpavus profectus, and the extant eulipotyphlan “insec- 
tivores” Solenodon paradoxus (Solenodontidae), Erinaceus europaeus (Erinaceidae), and Sorex araneus (Soricidae). We deleted all 
other non-chiropteran terminals. 

To this reduced, chiropteran-focused matrix, we added Vielasia sigei and a further 15 early Eocene bats Ageina tobieni, Archae- 
onycteris trigonodon, Ar. brailloni, ?Ar. storchi, Cambaya complexus, Eppsinycteris anglica, Hassianycteris messelensis, Honrovits 
joeli, H. tsuwape, Icaronycteris index,?!. sigei, Microchiropteryx folieae, Onychonycteris finneyi, Palaeochiropteryx tupaiodon, Proto- 
nycteris gunnelli; for family assignment see list of Comparative material, above). We did not include fossil species whose taxonomy is 
under review (e.g.,?/. menui which appears to represent a mixed species sample; SJH in prep.). Our final matrix therefore comprised 
22 chiropteran ingroup terminals (six extant bats, and 16 fossil bats representing all five currently recognised early Eocene families”), 
and eight non-chiropteran outgroup terminals (three extant, five fossil). Percentage completeness of our terminals for the cranioden- 
tal matrix is given in Data S1B. In the final 2665 craniodental character matrix, 925 characters (34.7%) were parsimony informative, 
and a further 433 characters (16.2%) were variable but non-parsimony informative. Following O’Leary et al.,°' all characters were 
specified as unordered. 

Maximum parsimony analyses of morphology matrix 

We analysed the full craniodental matrix using maximum parsimony (MP), as implemented in TNT 1.5,''® using a two-stage search 
strategy: an initial “New Technology” search with Sectorial Search, Ratchet, Drift, and Tree Fusing (with default settings for each of 
these) that was run until the same minimum tree length was found 100 times, followed by a “traditional” search using tree bisection 
reconnection within the trees saved from the first stage, saving a maximum of 100,000 trees. Multiple most parsimonious trees were 
summarised using strict consensus, and support values were calculated for clades in the consensus topology using 2000 bootstrap 
replicates, with bootstrap values reported as absolute frequencies. 

This analysis recovered two most parsimonious trees (length = 3165 steps), the strict consensus of which is shown in Figure S2A. 
Relationships among extant bats found in this analysis are partially incongruent with current molecular evidence: specifically, the 
pteropodid Pteropus was placed as the first branch within Chiroptera (rather than in a clade with the rhinolophoid Rhinopoma), sepa- 
rate from a clade comprising the remaining, laryngeally-echolocating extant bats (“‘microbats”) — a basal position for pteropodids is 
characteristic of most unconstrained morphology-only analyses of chiropteran phylogeny.°***° In this analysis, the fossil Eocene bats 
form a paraphyletic assemblage that is more closely related to the clade of laryngeally-echolocating extant bats than to Pteropus. 
Notably, Vielasia and Cambaya (known from a single 54 Ma dentary from the Cambay Formation at Vastan mine, India, that has 
been tentatively referred to Hassianycteridae**'*°) form a weakly supported (bootstrap <50%) clade, and this in turn is weakly sup- 
ported (bootstrap < 50%) as sister to the extant echolocating clade. Monophyly of the extant echolocating clade, to the exclusion of 
Pteropus and all Eocene fossil bats, is strongly supported (bootstrap = 97%). Besides the basal position of Pteropus, the overall to- 
pology in this analysis is broadly similar to our dated total evidence analysis (Figure 2). 

The fact that maximum parsimony analysis of the full craniodental matrix results in a topology that conflicts with current molecular 
evidence in placing Pteropus as the sister to all other bats (extant and fossil), rather than in a clade with our representative rhinolo- 
phoid, Rhinopoma (e.g., Shi and Rabosky®* and Amador et al.°“), raises questions about the ability of this matrix to accurately resolve 
bat phylogeny. To further investigate this, we filtered the matrix to retain only those characters that are fully compatible (= show no 
homoplasy) with the current molecular consensus view of bat phylogeny, as follows. We deleted all fossil taxa from the matrix, and 
then analysed the extant-only matrix using MP in TNT, using the same two-stage search strategy as before, and applying a “molec- 
ular scaffold” that enforced the monophyly of the following clades, based on recent molecular studies (e.g., Meredith et al.,? Álvarez- 
Carretero et al., and Upham et al. '7°): Eulipotyphla (our extant outgroup taxa Sorex, Solenodon and Erinaceus), crown Chiroptera 
(our six extant bat terminals), Yinpterochiroptera (Pteropus and Rhinopoma), Yangochiroptera (Saccopteryx, Nycteris, Myotis, and 
Pteronotus), and Emballonuroidea (Saccopteryx and Nycteris). This procedure resulted in a single most parsimonious tree. We 
then used the R package TreeSearch'”* to calculate the minimum possible length of each craniodental character on any tree (using 
the MinimumLength function), and its length on the single most parsimonious tree that resulted from the molecular scaffold analysis 
(using the CharacterLength function); if the latter was greater than the former (indicating some homoplasy relative to the molecular 
consensus), then that character was removed from the craniodental matrix. This procedure resulted in the removal of 519 characters, 
leaving 2146 characters, of which 467 were parsimony informative, and 422 were variable but non-parsimony informative. The fossil 
terminals were then re-added to this homoplasy-filtered matrix, which was analysed using MP in TNT with the same two-stage search 
strategy used for the full matrix. 

MP analysis of the homoplasy-filtered craniodental matrix resulted in three most parsimonious trees (length = 1598 steps). The 
strict consensus of these (Figure S2B) recovers relationships among extant bats that are fully congruent with the current molecular 
consensus, including monophyly of Yangochiroptera (bootstrap = 89%) and Yinpterochiroptera (bootstrap = 50%), although only the 
former is strongly supported. Similarly to MP analysis of the full craniodental matrix (Figure S2A), this analysis places all the Eocene 
fossil bats, including Vielasia, outside the crown clade, but bootstrap support for the crown clade is low (bootstrap <50%). The fossil 
bats are recovered in a single clade that is moderately strongly supported (bootstrap 73%), rather than as a paraphyletic assemblage 
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as found in the MP analysis of the full craniodental matrix (Figure S2A); within this fossil clade, Vielasia and Cambaya form a weakly 
supported (boostrap < 50%) clade that is sister to another weakly supported (bootstrap < 50%) clade comprising the remaining fos- 
sil taxa. 

Although this analysis similarly supports the hypothesis that Vielasia is not a crown clade bat, it is incongruent with previous 
studies**° (and with our maximum parsimony of the full craniodental dataset [Figure S2A] and our Bayesian dated analysis of our 
total evidence data set [Figure 2]) in placing all Eocene fossil bats in a clade, with onychonycterids (species of Onychonycteris, Ep- 
psinycteris, Ageina and Honrovits) in a nested position; these other analyses instead recover Eocene fossil bats as paraphyletic with 
respect to the crown clade, and place onychonycterids (including the well-preserved Onychonycteris, which shows a range of fea- 
tures that have been generally interpreted as plesiomorphic relative to other known fossil bats) in a basal position.° In addition, this 
kind of homoplasy-filtering approach represents a form of “maximum compatibility” analysis, which is a method of phylogenetic 
analysis that has had a long history of criticism (e.g., Farris and Kluge, ° Farris,'’° and de Pinna'’’'’°"'”), and characters that 
are homoplastic within one part of a phylogeny (e.g., within the crown clade) may still be useful for resolving relationships elsewhere 
in the tree, '’° and so arguably should not simply be discarded. For these reasons, we have used the full, unfiltered craniodental matrix 
for our total evidence analyses (see below). Nevertheless, we consider the results of this homoplasy-filtered craniodental analysis to 
provide additional evidence that Vielasia is not a crown clade bat. 

Templeton tests of alternative topologies 

To test whether our full and homoplasy-filtered craniodental matrices are compatible with Vielasia being a crown chiropteran, we 
carried out constrained MP analyses of both matrices with Vielasia enforced to fall within the crown-clade. As already noted, MP anal- 
ysis of the full craniodental dataset placed extant laryngeally-echolocating bats in a single clade to the exclusion of Pteropus, which 
was sister to all other chiropterans (Figure S2A); within the extant echolocating clade, Myotis and Pteronotus formed one subclade, 
and Nycteris, Rhinopoma and Saccopteryx the other. For the constrained analyses of the full craniodental dataset, we therefore en- 
forced Vielasia to either form a clade with Myotis and Pteronotus to the exclusion of all other taxa, or a clade with Nycteris, Rhinopoma 
and Saccopteryx to the exclusion of all other taxa. MP analysis of the homoplasy-filtered craniodental dataset (Figure S2B) recovered 
monophyly of Yangochiroptera and Yinpterochiroptera, and so for the constrained analyses of this dataset, we enforced Vielasia to 
either form a clade with the extant yangochiropteran taxa (Saccopteryx, Nycteris, Myotis, and Pteronotus), or the yinpterochiropteran 
taxa (Pteropus and Rhinopoma), to the exclusion of all other taxa. The constrained maximum parsimony analyses used the same 
search settings in TNT as the unconstrained analyses. 

Templeton tests '’° comparing the most parsimonious topologies from the unconstrained and constrained analyses were then im- 
plemented in PAUP* 4.0a169,'° the resultant p-values converted from two-tailed to one-tailed (following Goldman et al.'°°) by 
dividing by two, and mean p-values were calculated from the MPTs for each constrained analysis. In both cases, a constrained po- 
sition for Vielasia within the crown clade (or, in the case of the full craniodental analysis, within the extant echolocating clade) was 
rejected, with p-values for the Templeton tests <0.05 (Data S1E). Thus, a position for Vielasia within crown Chiroptera is strongly re- 
jected by our full and homoplasy-filtered craniodental matrices. 

Total evidence (TE) matrix 

Numerous studies support the use of a total evidence (combined morphological and molecular data) approach over morphology-only 
and molecular scaffold approaches for phylogenetic analyses that combine fossil and extant taxa (e.g., de Queiroz and Gatesy, '°' 
Manos etal.,'° Davalos et al.,'°° and Darlim et al. 184). We therefore created a total evidence matrix by combining our full craniodental 
matrix (see above) with molecular data for our six extant bat terminals and three extant non-bat outgroup terminals (Solenodon, Eri- 
naceus, and Sorex) taken from the alignment of O’Leary et al.,°' which comprises 36860 base pairs of DNA sequence data from 27 
nuclear loci, originally from Meredith et al.? All our fossil terminals were scored as unknown (?) for the molecular data. 

Undated Bayesian analysis of TE matrix 

Undated Bayesian analysis of the total evidence analysis was carried out in MrBayes v. 3.2.7a running on JASMIN, the UK’s collab- 
orative data analysis environment (hitps://jasmin.ac.uk'®°). We assigned a single Lewis '°° Mk model to the craniodental characters; 
given the large proportion of variable but non-parsimony informative characters (see above), we corrected for ascertainment bias by 
using the Mkv variant of the model, which assumes that only variable characters were scored (and so all constant characters were 
ignored in each analysis). An eight category lognormal distribution was used to model rate heterogeneity between characters. '°’ For 
the DNA sequence data, an optimal partitioning scheme and set of models was identified using PartitionFinder 2.2,''” with the 
sequence data initially partitioned by gene and (if protein-coding) codon position. The “greedy” algorithm and the Bayesian Informa- 
tion Criterion were used for model selection, with the choice of models restricted to those implemented by MrBayes, and with models 
that include both a gamma distribution and a proportion of invariant sites not tested'®° (Data S1C). 

The undated total evidence analysis comprised two independent runs of four chains with default heating parameters, run for 10 
million generations, sampling trees and other parameters every 5000 generations. To reduce the time needed to reach stationarity 
and convergence between chains, the following uncontroversial clades were enforced as monophyletic: Chiroptera (= all of our 
extant and fossil bat terminals), Eulipotyphla (= Solenodon + Erinaceus + Sorex), and Laurasiatheria (= all our terminals except Mae- 
lestes, Ukhaatherium, and Zalambdalestes). Tracer was used to identify an appropriate burn-in fraction for each analysis, and the 
post-burnin trees were summarised in MrBayes using 50% majority rule consensus, with Bayesian posterior probabilities as support 
values. 

The undated total evidence analysis (Figure S2C) provides strong support (BPP = 0.95) for monophyly of crown clade bats to the 
exclusion of all Eocene fossil bats included in the matrix, including Vielasia. The topology resembles that of the homoplasy-filtered 
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craniodental matrix (Figure S2B) in placing all Eocene fossil bats in a relatively strongly supported (BPP = 0.86) clade, within which 
Vielasia forms a strongly supported (BPP = 0.97) clade with Cambaya, and Vielasia + Cambaya is sister to the remaining fossil bats, 
with Onychonycteridae in a nested position. We consider this topology less plausible than that found in the dated total evidence anal- 
ysis (Figure 2) for reasons already discussed (see above), but it again supports the hypothesis that Vielasia is not a crown clade bat. 
Dated analysis of TE matrix 

Based on simulated data, Mongiardino Koch et al.°* found that including temporal information in Bayesian analyses using a fossilised 
birth-death (FBD) model led to increased phylogenetic accuracy relative to equivalent undated approaches. For these reasons, our 
preferred phylogenetic analysis is a dated Bayesian analysis of the total evidence matrix, and we focus on this here. 

For this analysis, we used the same substitution models for the morphological and molecular partitions as in the undated analysis 
(see above). We implemented the Fossilised Birth Death model with a combination of tip and node calibrations, following the general 
approach of Kealy and Beck, ‘°° Beck and Taglioretti, ‘°° and Beck et al.'°’ An Independent Gamma Rates (IGR) clock model’? was 
specified, with two separate clock models: one for the morphological characters, and one for the combined DNA sequence data. The 
prior on the clock rate was estimated using the custom R script of Gunnell et al.'°° A vague prior — uniform (0.001,200) - was specified 
for the variance of the IGR clock model, following Matzke and Wright. '°“ Diversity sampling was assumed, with a sample probability 
of 0.00418, based on the inclusion of six extant bat species out of a currently recognised total of 1437.'°° Because our taxon set does 
not include any fossil terminals younger than 47.4 Ma, we specified a shift in sampling at this time. '°° Priors on speciation, extinction 
and fossilisation were left as the MrBayes defaults. 

Extant terminals were assigned an age of 0 Ma, whilst fossil terminals were assigned a uniform age prior reflecting the stratigraphic 
uncertainty in their ages, following the recommendations of Püschel et al.'°’ (Data S1D). We enforced monophyly of Chiroptera, Eu- 
lipotyphla and Laurasiatheria (as in the undated analyses) and specified internal node calibrations for these three clades based on 
current fossil evidence, as follows. Chiroptera was assigned a minimum bound of 52.201 Ma (very slightly older than the minimum 
age of the oldest bat terminal in our matrix, Eppsinycteris anglica; Data S1D) and a maximum bound of 66 Ma. Eulipotyphla was as- 
signed a minimum bound of 62.5 Ma (the minimum age of Torrejonian To2, see Beck and Lee ‘°° and electronic suppl. mat. of Flynn 
et al.'°°) and a maximum bound of 83.6 Ma (the maximum age of the Campanian).”°° Laurasiatheria was assigned a minimum bound 
of 65.689 Ma (very slightly older than the minimum age of the oldest laurasiatherian terminal in our matrix, Protungulatum donnae) and 
a maximum bound of 113 Ma (the Aptian-Albian boundary”°"), which corresponds to the age of the oldest known eutherian with only 
four premolars (Sasayamamylos kawaii-°'), on the assumption that the origin of Placentalia (the last common ancestor of which is 
reconstructed as having four premolars®'), and hence Laurasiatheria, is highly unlikely to predate this. All three internal node calibra- 
tions were specified as offset exponential distributions,*°° such that there was a 5% prior probability for the calibrated node being 
older than the maximum bound. Preliminary analyses with the age of the root node unconstrained resulted in implausibly ancient pos- 
terior root ages, and so we used a truncated normal prior on the root age with a mean age (77.95 Ma) that was slightly older than the 
mean age (77.85 Ma) of the oldest terminals in our matrix (the outgroups Maelestes gobiensis, Ukhaatherium nessovi, and Zalamb- 
dalestes lechei), and a standard deviation of 1.0 Ma, following Sallam and Seiffert.2°° 

The dated total evidence analysis comprised two independent runs of four chains, with default heating parameters, run for 40 
million generations, and sampling trees and other parameters every 5000 generations. Stationarity and convergence between 
runs were determined using Tracer and used to identify an appropriate burn-in period; all post-burn-in trees were then combined 
using 50% majority rule consensus, with Bayesian posterior probabilities as support values. 

In the dated total evidence analysis (Figure 2), monophyly of crown-clade bats to the exclusion of all our fossil Eocene bats is 
strongly supported, with a Bayesian posterior probability (BPP) of 1.00. Relationships within the crown-clade are fully congruent 
with previous molecular studies, with very strong support for monophyly of Yangochiroptera (BPP = 1.00) and Yinpterochiroptera 
(BPP = 1.00). Vielasia sigei was recovered as the sister-taxon of crown-clade bats in this analysis, although with only moderate sup- 
port (BPP = 0.63); the presence of a non-zero branch length leading to V. sigei indicated that it is not plausibly ancestral to the crown 
clade. The sister-taxon to Vielasia + the crown-bat clade was Cambaya, and monophyly of Cambaya + Vielasia + crown clade bats 
received strong support (BPP = 0.95). 

This analysis also recovered a fairly strongly supported (BPP = 0.85) clade comprised of onychonycterids (species of Onychonyc- 
teris, Eppsinycteris, Ageina and Honrovits) and another containing the other four extinct Early Eocene bat families in our sample 
(icaronycterids + archaeonycterids + hassianycterids + palaeochiropterygids), although support for this latter clade is very weak 
(BPP = 0.51). Onychonycteridae was the first family to diverge, which is congruent with previous suggestions. ®™?®31:67 Within Ony- 
chonycteridae, relationships were fully resolved and consistent with previous studies.2°*' Archaeonycterids (Archaeonycteris and 
Protonycteris spp.) form a poorly supported (BPP = 0.51) clade, with the exception of A. trigonodon which is moderately strongly 
supported as sister to icaronycterids + hassianycterids + palaeochiropterygids (BPP = 0.82). Within this latter group, the relationship 
of Icaronycteris spp. to each other and to a clade comprising the hassianycterids + palaeochiropterygids were unresolved. The 
possible paraphyly of Archaeonycteris and Icaronycteris has been suggested previously.7°"*' 

Maximum parsimony optimisation of the craniodental matrix on the dated total evidence phylogeny (our preferred topology; Fig- 
ure 2) in PAUP 4.0a'° results in 49 characters being identified as unambiguous (i.e., under both accelerated and delayed transfor- 
mation) synapomorphies for crown Chiroptera to the exclusion of Vielasia; a full list of synapomorphies for all nodes in the dated total 
evidence phylogeny is given in Data S2. 
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Bayes factor tests of alternative topologies 

To test whether our total evidence matrix is compatible with Vielasia being a crown chiropteran, we carried out constrained versions 
of the undated and dated Bayesian analyses with Vielasia enforced to fall within the crown-clade; we then compared the fit of these 
constrained analyses to the unconstrained analyses using Bayes factors.7°* °°° Both undated (Figure S2C) and dated (Figure 2) anal- 
ysis of the total evidence dataset recovered monophyly of Yangochiroptera and Yinpterochiroptera, and so for the constrained ver- 
sions of these analyses, we enforced Vielasia to either form a clade with the extant yangochiropteran taxa in the matrix (Saccopteryx, 
Nycteris, Myotis, and Pteronotus), or the yinpterochiropteran taxa (Pteropus and Rhinopoma), to the exclusion of all other taxa. We 
calculated harmonic means for the unconstrained and constrained topologies using stepping stone analysis*” calculated using 10 
million generations, with otherwise default settings. Bayes factors were then calculated as 2x the difference in log likelinood between 
the constrained and unconstrained topologies.*°° °°” For all constrained topologies, Bayes factors (B10) for the stepping stone an- 
alyses were >150 (Data S1E), representing very strong evidence’ °°’ against a position for Vielasia within crown Chiroptera, and in 
favour of its unconstrained position as a stem chiropteran. 

Summary of phylogenetic analyses 

In summary, the results of the phylogenetic analyses presented here collectively provide very strong evidence that Vielasia is not a 
crown clade bat (see Figures 2 and S2A-S2C; Data S1E). In turn, they support the hypothesis that the advanced echolocation abilities 
evidently present in Vielasia based on its cranial anatomy originated outside the clade that encompasses all modern laryngeally- 
echolocating bats. However, our preferred topology (Figure 2), from the dated total evidence analysis, places Vielasia sister to the 
crown-clade, and so suggests that the echolocation abilities of Vielasia may be homologous with crown bats, and thus that the laryn- 
geal echolocation characteristic of modern non-pteropodid bats evolved once, prior to the diversification of the crown clade. 


LINEAR AND GEOMETRIC MORPHOMETRIC ANALYSES 

Phylogenetic framework 

In order to investigate the impact of phylogenetic non-independence on our morphometric analyses, we used phylogenetic compar- 
ative methods and three different dated phylogenies of bats: (1) the total evidence tip-dating phylogeny produced as part of this study 
(e.g., Figure 2), (2) the molecular node-dating phylogeny of Shi and Rabosky,®° and (3) the molecular node-dating phylogeny of Ama- 
dor et al.°* The total evidence tip-dating phylogeny obtained here only allows for taking into account variation in branch lengths, 
which are proportional to time: Vielasia sigei, as a tip, and the node where Vielasia diverges from the crown bat clade, are recon- 
structed with a range of ages based on the post-burnin trees output by MrBayes. The phylogenies of Shi and Rabosky*? and of Ama- 
dor et al.°* allow for taking into account variation in both topology and branch lengths (= divergence times). Combining these three 
different phylogenetic datasets, we could therefore take into account variation in both topology and estimated divergence times (the 
only constant point being the stem position of Vielasia sigei), and this served to test the robustness of our statistical results. 

Vielasia is present in our total evidence tip-dating phylogeny, as the sister taxon of crown Chiroptera, but needed to be added in this 
position to the Shi and Rabosky®* and Amador et al.°* molecular phylogenies a posteriori. This is only possible if the age of the 
Vielasia + crown Chiroptera node (as estimated in our total evidence tip-dating phylogeny) predates the age of crown Chiroptera 
in the two molecular phylogenies. We sampled 1001 post-burnin trees output from our total evidence tip-dating analysis and filtered 
them to retain only those in which the age of Vielasia + crown Chiroptera predates the age of crown Chiroptera estimated by Shi and 
Rabosky®° and Amador et al.°*; to do this, we used the packages treeio,”°° ape,”°° phytools,?'° and ULT,7'' which resulted in 104 
trees for the Shi and Rabosky®® phylogenetic dataset, and two trees for the Amador et al. dataset.°* 

However, this approach resulted in only a small amount of variation in branch lengths being retained (particularly in the case of 
Amador et al.,°* where only two trees can be used), and, perhaps more importantly, it resulted in a potentially unrealistically short 
branch between Vielasia + crown Chiroptera and crown Chiroptera, which would reduce the probability that the similarities between 
Vielasia and crown chiropterans are convergent: the branch separating Vielasia + crown Chiroptera from crown Chiroptera is so short 
that very little evolutionary change will be estimated to have occurred along it, and so similarities between Vielasia and crown chi- 
ropterans will be likely to be interpreted as due to common ancestry. Therefore, we also rescaled the Shi and Rabosky®’ and Amador 
et al. molecular phylogenies so that the age of crown Chiroptera matched the age of this node in the 1001 post-burnin trees ob- 
tained from our total evidence tip-dating analysis. This allows for a greater range of variation in branch lengths to be taken into ac- 
count. To do this, we used the R packages scales” '* and geiger’'* to linearly rescale the ages of the Shi and Rabosky®* and Amador 
et al. molecular phylogenies so that the age of crown Chiroptera matched the age of this node in each of the 1001 post-burnin trees 
from our total evidence tip-dating analysis, with the ages of more nested nodes adjusted by the same (linear) scaling factor. 

We note here that the phylogeny of Amador et al. includes a node that is younger than one of its two descendants; this is an arte- 
fact produced by the software used by the authors of that article (N.P. Giannini, pers. comm.). To avoid problems with downstream 
analyses, we changed the age of this node so that the length of the branch leading to its “older” descendant became half the minimal 
length of all other internal branches (i.e., half of 0.0648 Ma); this represents a compromise between keeping a negative value and 
replacing it with an arbitrary positive value (i.e., the branch length is now positive instead of negative, but still half the length or 
less than that of any other internal branch). In addition, we note that most of the extant tips in the phylogeny of Amador et al.** do 
not have an exact age of 0 Ma (out of 812 taxa, only four have an age of 0 Ma, 30 have an age of 0.24 Ma, and 778 have an age 
of 0.35 Ma). However, we did not correct these very minor anomalies, as they would have a minimal impact on our downstream 
analyses. 
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As a result of these procedures, we were left with 1105 trees corresponding to the Shi and Rabosky® phylogeny — 104 with original 
ages that are already compatible with the post-burnin trees from our total evidence tip-dating analysis, and 1001 with ages rescaled 
to match the post-burnin trees -, and 1003 trees corresponding to the Amador et al.°* phylogeny — two with original ages that are 
already compatible with the post-burnin trees, and 1001 with ages rescaled to match the post-burnin trees. In both cases, two of 
the rescaled trees had a null branch length between the Vielasia + crown Chiroptera node and the crown Chiroptera node; these 
two trees were removed, leaving 1103 trees corresponding to the Shi and Rabosky®* phylogeny, and 1001 corresponding to the 
Amador et al. phylogeny.®* 

Bony labyrinth morphometric analysis 

Using several morphometric features of the bony labyrinth, Davies et al.'°°° were able to discriminate between bats and non-bat 
placental mammals, as well as between echolocating and non-echolocating bats. Of these features, we focused on two we were 
able to measure on the bony labyrinth of Vielasia sigei and that may give an indication of the echolocating ability of this species. 
The first is the length of the cochlear basilar membrane, the allometric relationship of which with body mass led Davies et al.'° to 
note that, in comparison to other placentals, echolocating bats have a longer cochlear basilar membrane than expected given their 
body mass. Secondly, Davies et al.°° calculated the size of the semicircular canals (estimated as the average of a canal’s height and 
width) and the size of the cochlea (a one dimensional value calculated following Spoor et al.°’ as the average of the diameters of the 
first and second cochlear turn and of the “slant height” of the cochlea) and found a difference between echolocating bats and other 
mammals (including non-echolocating bats) when plotting semicircular canal size against cochlea size (echolocating bats having 
larger cochlea/smaller canals). 

For both basilar membrane length vs body mass and semicircular canal size vs cochlear size, we tested whether we could discrim- 
inate between bats and non-bat placentals, and between different echolocation types in bats (non-laryngeal echolocating, oral-emit- 
ting, and nasal-emitting). We included cetaceans in both datasets because they represent a mammalian clade that includes species 
that have evolved to hear very high frequencies (odontocetes), similarly to bats, but also species that have evolved to hear very low 
frequencies (mysticetes). These two additional classes may allow us to identify a functional signal in the traits tested and therefore 
increase the validity of our inferences. After controlling for potential effects of these ecological categories on the relationships we 
tested, our final goal was to include Vielasia and test whether it could be confidently assigned to a particular class (i.e., non-LE 
bats, LE-bats, HF-hearing cetaceans, LF-hearing cetaceans, or other placental mammals). 

For the cochlear basilar membrane length vs body mass relationship, we simply added measurements for Vielasia to the dataset of 
Davies et al. (see table S2),'° with therefore 82 species included (44 extant bats, 37 extant mammals, and Vielasia sigei; Data S1G). 
The dataset of Davies et al. (see table S1)°° included measurements of cochlear and semicircular canal sizes for chiropterans only; 
thus, we had to obtain the non-chiropteran data they used directly from various references,°*"°*"''* to which we added additional 
measurements for cetaceans.''”''' The resulting dataset therefore includes 124 species, with 52 extant bats, 71 extant mammals, 
and Vielasia sigei (Data S1M). The phylogeny used by Davies et al.°° for their comparative analyses does not include cetaceans, and 
so we instead used the dated mammalian phylogeny of Alvarez-Carretero et al.° as the starting tree for our own analyses. However, 
we replaced the bat clade in the Alvarez-Carretero et al. phylogeny with the bat phylogenies derived from Shi and Rabosky,®° Ama- 
dor et al.,°* and the current study (see “phylogenetic framework” above). This requires that the age of Vielasia + crown Chiroptera is 
younger than that of the node uniting Chiroptera with its sister taxon, which in the Alvarez-Carretero et al.° phylogeny is Artiodactyla + 
Perissodactyla + Carnivora + Pholidota. This was the case for 74 of the 104 trees derived from Shi and Rabosky** without rescaling, 
but was not the case for either of the two trees derived from Amador et al.°* without rescaling (see “phylogenetic framework” above 
for details); for the rescaled phylogenies, this was the case for 971 of the 1001 trees derived from Shi and Rabosky,°° and for 974 of 
the 1001 trees derived from Amador et al.°* For each dataset (basilar membrane length vs body mass; semicircular canal size vs 
cochlear size), we then pruned the relevant phylogenies to retain only those species for which we had measurements. 

Basilar membrane length vs body mass - Regarding the allometric relationship of basilar membrane length to body mass, we first 
simply plotted all points, distinguishing them according to the following groups: high-frequency hearing cetaceans, low-frequency 
hearing cetaceans, non-laryngeally echolocating bats (pteropodids), laryngeally echolocatng bats (nasal- and oral-emitters), and 
other placentals. 

We then performed phylogenetic generalized least squares (PGLS) regressions of (natural logarithm of) basilar membrane length 
against body mass, repeating the analysis for all phylogenies we created. To do so, we used functions and methodology from the 
packages phylolm?'* and sensiPhy*'® to implement PGLS regression of log basilar membrane length against log body mass, taking 
into account phylogenies with variation in their relationships (within Chiroptera) and in their branch lengths. The function phylolm from 
the eponym package’ “ allows for performing PGLS regressions with several different possible evolutionary models (Brownian mo- 
tion, Brownian motion with trend, Ornstein-Uhlenbeck with a fixed and estimated root value or with a stationary ancestral distribution, 
Pagel’s lambda, kappa, and delta, and the early burst model). The function tree_phylm from the sensiPhy package enables the PGLS 
regression to be repeated across multiple phylogenies, and outputs tables to compare averaged parameters between the tested 
models. We wrote a function named tree_phylm_full_results'*° which is a slightly modified version of the tree_phylm function (of 
the sensiPhy package). This custom function can be easily further modified to output more information than given by the standard 
tree_phylm function, such as regression residuals, fitted values etc. This custom function returns two important elements that are 
also returned by the standard tree_phylm function, namely: the estimated parameters for each PGLS regression (the $sensi.esti- 
mates element) and the averages of the parameters for each model (the $all.stats element). In tree_phylm_full_results however, 
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the $sensi.estimates element as well as other outputs from the regressions can be returned for all tested models, not just for the best- 
fitting model, the choice being up to the user (by turning the option output.all.models to TRUE or FALSE). 

Finally, we tested whether there were significant differences in the PGLS regression parameters between selected comparisons 
(bat versus non-bat placentals, laryngeally echolocating bats versus other placentals, high-frequency hearing mammals versus 
others, high- and low-frequency hearing mammals versus others). To perform these comparisons, an additional wrapper on the 
tree_phylm function, compare.tree_phylm'° was written. This function takes the same input variables as the tree_phylm_full_results 
function, plus information regarding potential subsetting of the dataset. If a subset is specified, the dataset is divided in two groups 
(i.e., both the phylogenies and the morphometric data), PGLS is performed for both (using tree_phylm_full_results function), then the 
compare.tree_phylm function provides the results from a t-test contrasting average slopes and intercepts for the PGLS with best- 
fitting evolutionary models for each group, and finally it outputs average regression line parameters to be able to simply add them 
to an existing biplot. 

Overall, there are few differences between our comparisons (Data S1H), that aim to compare selected subsets of our dataset (LE 
bats, all bats, LE bats plus HF-hearing cetaceans, LE bats plus both HF-hearing and LF-hearing cetaceans) both to the remaining 
species (i.e., non LE bat placentals, non-bat placentals etc.) and to the whole dataset. There is no significant difference in regression 
slopes (p > 0.4 for every comparison), and only some significant differences in regression intercepts. These differences generally 
contrast echolocating mammals and other species: LE bats seem to differ from other placentals, although this relationship is not sig- 
nificant (p = 0.06), bats are significantly different from non-bat placentals (p = 0.02), and placentals with extreme hearing abilities (i.e., 
HF-hearing and LF-hearing mammals) in general differ markedly from other placentals (p = 0.009). 

No selected subset differs significantly from the overall placental trend (i.e., the regression for the whole placental dataset), and so 
we only added the placental regression line (i.e., Pagel’s delta with an average value of 2.47) to the primary drawn biplot in the main 
text (Figure 3B); other configurations can be found in Figure S3A. Regarding the position of Vielasia, it is located in the middle of the 
morphospace comprising laryngeally echolocating bats, which is far from that of non-echolocating bats. This suggests echolocating 
abilities in Vielasia, although echolocating bats do not differ enough from other mammals to demonstrate this statistically. Morpho- 
spaces of both echolocation call types (i.e., emission through the nose or the mouth) overlap considerably, and Vielasia falls within 
this area of overlap; thus, this particular analysis is unable to determine the type of echolocation emission of Vielasia. 

For our analysis of the relationship between semicircular canal size vs cochlear size, we first performed a simple, graphical com- 
parison of semicircular canal size (radius of their curvature, R) and cochlea size (as defined by Spoor et al., but see also Davies 
et al.°°) for bat species as done by Davies et al.°° We however diverged from the original representation of Davies et al. (see figure 4)°° 
by instead performing boxplots (to better represent whole bat trend) adding values ranges (for LE yangochiropterans, LE yinptero- 
chiropterans, and non-LE yinpterochiropterans) to compare with the values of Vielasia sigei (Figure 3C). Values were not corrected for 
relative size or phylogeny. 

Second, our goal was to replicate Davies et al. regressions of each semicircular canal size against cochlea size (e.g., see figure 
5 in Davies et al.°°). Instead of analyzing each regression separately, we aimed to encompass all three regressions and provide a 
multivariate depiction of the relationship between semicircular canal size and cochlea size. We could not directly perform multivariate 
analyses, as the sensiPhy package*'° cannot deal with multivariate data, and as there is currently no satisfactory method for multi- 
variate PGLS - there are some methods for highly correlated variables that are landmark coordinates, but as yet not for putatively 
independent ones (e.g., Adams*'° and Clavel and Morlon?'’). Therefore, we performed a PGLS regression for each semicircular ca- 
nal separately (i.e., natural logarithms of semicircular canal radius against cochlea size), and PGLS residuals were then used as the 
basis for further interpretation. It should be noted that these residuals are not phylogeny-corrected residuals; they are simply more 
appropriate to use than residuals derived from OLS regressions.’ Moreover, subsequent tests on regression residuals are not rec- 
ommended since they can give erroneous results.?!™??? Hence, we took an appropriately cautious approach, applying descriptive 
treatments: we performed discriminant analyses taking into account phylogeny (‘phylogenetic flexible discriminant analysis’ or 
pFDA") or not (traditional linear discriminant analysis, LDA, using the package MASS°°') to examine first whether the link between 
canal and cochlea size proposed by Davies et al.°° depends on phylogeny (see Data S11), and second in which group Vielasia sigei 
falls in each case (see Data S1J). To determine whether Vielasia could be assigned to a particular hearing ability group (same groups 
as the basilar membrane length vs body mass relationship), we specified the groups as the discriminant factor, Vielasia as the “test” 
group, and all other extant species as the “training” set of the discriminant analysis (as their membership of a particular group is 
already known). We repeated these analyses for all the phylogenies we created, and we reported average coordinates, reflecting 
the coordinate variation for each species, for each analysis (Figures 3D and S3B). 

Numerous R packages allow for performing PGLS regression; here, we chose to use the pgls function of the package caper.*”* The 
pFDA function (to perform the eponym analysis), as written by Motani and Schmitz, °? takes phylogeny into consideration by modi- 
fying the inputted phylogeny by transforming its branch lengths (and therefore its variance-covariance matrix) to account for phylo- 
genetic autocorrelation using Pagel’s lambda. Motani and Schmitz'*° recommend estimating the lambda parameter from the data, 
and to provide this estimate of lambda to the pFDA function. However, as well as values of Pagel’s lambda, the pgls function of the 
package caper estimates (using maximum likelihood) values for two other models: Pagel’s kappa and delta (see for instance the 
vignette of the package caper”). By default, Pagel’s kappa and delta are not taken into account in the pFDA function. However, 
since all three models are direct transformations of the tree, we slightly modified the pFDA function written by Motani and Schmitz‘? 
to also take these two parameters into account, applying the exact same transformation to the tree that the function pgls does. 


, sê? 
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Doing three regressions (i.e., each semicircular canal size vs cochlear size), the three values for each Pagel’s parameter are likely to 
differ, whereas the pFDA function can only handle a single value for each (i.e., a single way to transform the tree). A unique value for 
each parameter has therefore to be found that can be applied to all three regressions; we wrote and used the function conv.pgls to 
deal with this.'*° This function runs PGLS regressions of each variable, narrowing progressively the interval in which the values are 
estimated, which results in parameter estimates that converge asymptotically. First, the conv.pgls function performs PGLS regres- 
sions with the default bounds for each parameter; each of these first regressions yield 95% confidence intervals for each parameter. 
Then, the function retains the shared interval of all 95% intervals for each parameter, and repeats PGLS regressions, narrowing pro- 
gressively this interval using a user-defined thinning parameter (with a default value set to 10%) that is not necessarily symmetrical: 
the narrowing of this interval of a (provided) % value is proportional to the position of the parameter value (i.e., the closest to a bound 
is the parameter value, the less the interval is narrowed on that bound, and the more it is narrowed on the other bound, the sum of this 
narrowing being equal to the provided % value). At the end, the function outputs the “optimal” common value of the three parameters 
lambda, kappa, and delta, together with two tables comparing 1) AIC values for regressions with locally (i.e., for each regression) or 
globally (i.e., for all regressions) optimal parameters, and comparing 2) the parameters value and 95% confidence interval of each 
PGLS regression with the “globally optimal” parameter values. The function can also output the residuals for each PGLS with the 
“globally optimal” parameter values. We then wrote two functions to repeatedly perform pFDA and LDA (pfda_LSchmitz and rep_Ida 
respectively'*°) using the “globally optimal” Pagel’s parameter values (and the residuals arising from PGLS regressions using them) 
for all phylogenetic trees considered. These functions output, for each iteration (= phylogeny), a confusion matrix (a matrix of true and 
predicted attributions for the training dataset) and various results regarding the training and test datasets (i.e., predicted class, prob- 
ability for each class, discriminant axes scores). 

To help visualize the variation between all calculated pFDAs and LDAs (i.e., depending on phylogenies) for each species’ score on 
discriminant axes, we also plotted the morphospaces of each species’ points, using both all points of each species (‘full’ morpho- 
spaces) and the 95% points closest to the ‘full’ morphospace centroid (‘'95%’ morphospaces; i.e., points below the 95% quantile of 
Euclidean distance to the full morphospace centroid, the latter being calculated using the sf package**°). This procedure has been 
wrapped in a custom function named plots. variation, "°° taking into account data points (and morphospaces if already computed and 
saved), average positions, analyses names, axes information (number, percentage of contribution), ecological groups (and colors) if 
applicable, and options to return the computed morphospaces and to output the produced figures. These plots of ‘full’ and ‘95%’ 
morphospaces (Figure S3B) allow for more detailed comparison of the position of Vielasia relative to the morphospaces of the other 
taxa in these plots, and to visualize the full variation in discriminant axes scores. 

In both pFDA and LDA, the first discriminant axis accounts for ~75% of the between-group variance. There is also a clear link be- 
tween phylogeny and hearing ability type across placentals. Though there are errors in the classification (Data S11), there is a hier- 
archy along the first discriminant axis in the LDA (i.e., when phylogeny is not taken into account), with high-frequency hearing non-bat 
mammals (odontocetes), then echolocating bats plus low-frequency hearing non-bat mammals (mysticetes), and then non-echolo- 
cating bats and all other placentals (Figures 3D and S3B). Most of the classification errors are swaps between overlapping groups. In 
the LDA, Vielasia is classified as an echolocating bat (Data S1J; Figures 3D and S3B), but nasal-emitting and oral-emitting bats are not 
distinguishable along the first linear discriminant axis. Taking into account phylogeny by means of pFDA, there seems to be no 
ecological signal at all: all taxa are predicted as “non-echolocating and non-bat placentals” (Data S1J), and all categories overlap 
(Figure S3B). Therefore, Vielasia is also retrieved as a “non-echolocating and non-bat mammal” (Data S1J) and falls roughly in 
the middle of the plot of the two first discriminant axes (Figure S3B). Variation in the phylogeny does not alter these observations: 
there is more variation in species position in the pFDA than in the LDA (Figure S3B), but it remains slight. 

In summary, the echolocation call types/hearing abilities may depend on phylogeny, phylogenetic variation not much altering 
global trend; we hence preferred to provide average species’ position and groupings of the LDA in the main text (Figure 3D), detailed 
graphs for both analyses being available in the Figure S3B. According to these results, Vielasia would presumably have been able to 
echolocate because it is a bat, but non-LE bats still fall within the morphospace of “other” placentals; in bats, phylogeny does not 
exactly “predict” echolocation call type, and Vielasia’s position in the LDA may be independent from phylogeny as well. 
3D geometric morphometric analysis 
We included the composite 3D cranial and mandibular models of Vielasia sigei in the 3D geometric morphometric dataset of Arbour 
et al.,°° E®M3 capturing shape of bat crania (with 36 landmarks and 5 curves represented by equidistant semi-landmarks) and man- 
dibles (with 20 landmarks and 4 curves). Using the same protocol, we collected these data from Vielasia sigei with Landmark 3.6''* 
(Data S1N and S10). This allowed us to compare Vielasia with 202 extant bat species for the cranium and 191 species for the 
mandible, using the dietary and echolocation emission categories assembled by Arbour et al.°° (Data S1P). Processing of the cra- 
nium and mandible landmark data (i.e., inputting the data, flipping the data when the left or right landmark was missing, estimating 
missing data, performing the generalized Procrustes alignment, averaging Procrustes coordinates for each species, and averaging 
paired landmarks) was performed in R 4.1.1, '' following the approach of Arbour et al.,°° and was implemented as a custom R func- 
tion (Arbour_et_al_treatment; available at Maugoust et al. '*°); this function uses functions from the packages geomorph??4??S and 
Losi 

We then performed phylogenetic Principal Component Analyses (pPCA) on these aligned datasets using the package phytools. In 
order to examine variation in pPCA results implied by phylogenetic variation (in terms of branch lengths), pPCAs were performed for 
all 1103 trees corresponding to the Shi and Rabosky?? phylogeny (see “phylogenetic framework” above), and so there were a total of 
1103 pPCAs for the cranial and the mandibular datasets. To reproduce Arbour et al.’s°° approach as closely as possible (and so 
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ensure meaningful comparison between our results and theirs), we did not use the set of phylogenies corresponding to the later phy- 
logeny of Amador et al.,°* as it was not used by Arbour et al.°° 

Parallelisation of the pPCAs computation was done using the packages foreac and doParallel.22° A problem of repeated 
pPCAs with slightly varying inputs we encountered is that the outputted principal components can be inverted in their sign depending 
on the inputs; a custom function (align_rep_ppca) was therefore written to resolve this issue (see comments within this function for 
technical details'~°). 

After this procedure, average pPCA coordinates were calculated for each species on each phylogenetic Principal Component 
(pPC) axis over all 1103 iterations. As done by Arbour et al.,°° the proportion of variance explained by each pPC axis was calculated 
using the procD.Im function of the geomorph package; here this was done for each iteration and then averaged. We then plotted the 
average coordinates of the two first axes of the pPCAs for cranial and mandibular landmarks (Figure 4). 

We also plotted the variation in species’ position for all pPCA axes (axes 1-3 for the cranium and 1-4 for the mandible) with ‘full’ and 
‘95%’ morphospaces (see previous section) together with average coordinates (Figure S5). Regarding the cranium (Figure S5A), 
there is relatively little spatial variation within the ‘full’ and ‘95%’ morphospaces for each species, with the most variation on the sec- 
ond and third pPC axes. The most significant point is that, on a plot of pPC1 and pPC2 (upper row of Figure S5A), Vielasia plots well 
outside the morphospace of non-echolocating bats, and within the morphospace of oral-emitting echolocating bats, with a clear 
offset between Vielasia and non-echolocating bats on both axes. Regarding the mandible (Figure S5B), the morphopace for each 
species is once again quite small. In plots of pPC1 and pPC2 (upper row of Figure S5B), most notable is the fact that the morphospace 
of carnivorous bats overlaps those of nectarivorous and omnivorous bats. In plots of pPC3 and pPC4, each species has a slightly 
larger morphospace, especially for sanguivorous bats. However, Vielasia’s position within dietary morphospace remains similar: it 
is near the middle of each pPC axis, corresponding to the overlapping morphospaces of carnivorous, frugivorous, and insectivorous 
bats on pPC1 and pPC2, and at the edge of the insectivorous morphospace in pPC3 and pPC4. 

Underlying methodological implications 

The goal of our set of morphometric analyses described above was, beyond their original purpose, to test how the results were 
affected when phylogenetic uncertainty (in terms of phylogenetic relationships and of branch lengths) was taken into account. We 
first found that, after repeating PCAs (arising from geometric morphometrics) and DAs (of regression residuals) over all the phylog- 
enies we considered, exhaustive graphical results (i.e., showing the whole morphospace of each species) do not differ from simpler 
ones (i.e., showing only the average values of each species coordinates); this justifies our use of the simpler graphs in the main text 
(Figures 3B, 3D, and 4). Secondly, again using information from phylogenetic variation, we found little difference between echolocat- 
ing and non-echolocating placental mammals regarding the allometry of cochlear basilar membrane length (Figure S3A; Data S1H). 
Alternative allometries of cochlear basilar membrane length (Figure S3A; Data S1H) show that the distinction between laryngeally 
echolocating bats and other placental mammals found by Davies et al.'° may be less robust than first thought, but there would still 
be an offset of the allometric relationship between placentals with extreme hearing abilities and the other species. In all cases, Vielasia 
lies almost on the regression line of the group of interest. Finally, phylogeny is likely to play an important role in ecological discrim- 
ination in the relationship between semicircular canal size and cochlear size (Figure S3B; Data S11 and S1J), and should be consid- 
ered when inferring the ecology of extinct species. 


h22" 


QUANTIFICATION AND STATISTICAL ANALYSIS 


The statistical treatment of morphometric data was performed using R 4.1.1'*' and PAST 4.05''° (see key resources table). It can be 
re-performed using the source data (also available in Data S1G, S1K, S1M-S1P), with the R script and the entire results of its execu- 
tion (a RData file) being available at.'*° The R script contains original R functions written for our analyses (as detailed above) and ex- 
isting R functions (with minor modification) to perform phylogenetic flexible discriminant analysis ((FDA??; retrieved from the Github 
page of Lars Schmitz). The phylogenetic analyses were carried out in TNT 1.5''° (maximum parsimony analyses of the morphological 
matrices) and MrBayes 3.2.7a’’ (undated and dated Bayesian analysis of the total evidence matrix), with Templeton tests of alter- 
native topologies and mapping of synapomorphies in PAUP* 4.0a169'2°; the morphological and total evidence matrices as TNT- 
formatted and MrBayes-formatted files respectively, and the trees output by analyses of these matrices, are available at Maugoust 
et al.'°° DOls are listed in the key resources table. 
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